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Abstract 

We present consistently ordered calculations of the structure functions ^{x^Q 2 ) and 
J~l(x,Q 2 ), in different expansion schemes. After discussing the standard expansion in powers of 
ot s {Q 2 ) we consider a leading-order expansion in ln(l/x) and finally an expansion which is leading 
order in both ln(l/x) and a s (Q 2 ), and which is the only really correct expansion scheme. Ordering 
the calculation in a renormalization-scheme-consistent manner, there is no factorization scheme de- 
pendence, and the calculational method naturally includes to the "physical anomalous dimensions" 
of Catani. However, it imposes stronger constraints than just the use of these effective anomalous 
dimensions. A relationship between the small-x forms of the inputs J-2(x,Q 2 ) and J-l(x,Q 2 ) is 
predicted. Analysis of a wide range of data for J-2(x, Q 2 ) is performed, and a very good global fit 
obtained, particularly for data at small x. The fit allows a prediction for !Fl(x, Q 2 ) to be produced, 
which is smaller than those produced by the usual NLO-in-a s (Q 2 ) fits to T 2 ix,Q 2 ) and different 
in shape. 



July 1997 



1. Introduction. 

The recent measurements of J-2(x,Q 2 ) at HERA have provided data on a structure function 
at far lower values of x than any previous experiments, and show that there is a marked rise in 
F 2 (x,Q 2 ) at very small x down to rather low values of Q 2 [1][2]. Indeed, the rise persists for values 
of Q 2 as low as 1.5 GeV 2 . 

The qualitative result of a steep rise at small x was in conflict with standard methods used 
to fit the data were based on the solution of the Altarelli-Parisi evolution equation [3] at the two- 
loop level, using nonperturbative, flat (the Donnachie-Landshoff pomeron used to soft physics has 
behaviour x -0 08 [4], and we will take steep to mean any powerlike behaviour steeper than this) 
input parton distributions at starting scales of Q 2 ~ 4GeV 2 (e.g. [5]). This procedure results in 
an effectively steep behaviour at small x [6], but only after a long evolution length, and therefore 
at values of Q 2 3> 4GeV 2 . Thus, the data led to optimism amongst those advocating the use of 
the BFKL equation [7] , which provides the unintegrated gluon Green's function which includes the 
leading power of ln(l/x) for any power of (fixed) a s . It was solved analytically in the asymptotic 
limit x — > 0, or numerically for finite x, and predicted a powerlike behaviour of x~ 1_A for the gluon 
distribution function, where A = 41n2a s and a s = (3/ir)a B , i.e. A ~ 0.5 if a s ~ 0.2. This was 
assumed to lead to ^{xjQ 2 ) behaving like x~ x and was in rough qualitative agreement with the 
data, even if A was somewhat high. It could also be seen as some justification for choosing powerlike 
inputs (with A ~ 0.2 — 0.3) for the parton distributions (e.g. [8]), which could then enable a good 
fit to the data using the Altarelli-Parisi equation. 

However, it was also convincingly demonstrated that it was possible to generate the observed 
steep behaviour by being a little less conservative concerning the region in which perturbative 
evolution could be applied. Gliick, Reya and Vogt had in fact predicted a sharp rise in F 2 (x,Q 2 ) 
at small x, even for Q 2 ~ lGeV 2 , by using two-loop evolution from roughly valence-like parton 
distributions at a starting scale of Q 2 = 0.34GeV 2 [9]. Starting at a higher scale, Q 2 = lGeV 2 , Ball 
and Forte were able to fit the small-x data using their double asymptotic scaling (DAS) formula 
[10], which is a simple, but very accurate, approximation to the solution of the one-loop evolution 
equation with flat 1 inputs and which is valid in the region of small x (x & 0.01). They also showed 
that two-loop evolution with flat inputs starting at Q 2 « 2GeV 2 could fit the data available at 
that time very well [11]. 

This state of affairs left scope for argument about the real underlying physics describing the 
small-x behaviour of Ti{x,Q 2 ). Those who used Altarelli-Parisi evolution from small scales could 
be accused of working in regions where perturbation theory was questionable and, perhaps more 
importantly, of ignoring terms of higher order in a s (but also higher order in ln(l/x)), which 

1 Parton inputs behaving like a; -1 , as opposed to x~ 1_A , are referred to as fiat, i.e. the parton density 
rescaled by x is fiat. 
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seemed from the BFKL equation to have very important effects. Conversely, those who used 
the BFKL approach could be accused of ignoring all but the leading-ln(l/x) terms (and hence 
ignoring the large-x data) and also of working in a less well-defined theoretical framework than the 
renormalization group approach based on the factorization of collinear singularities [12]. Starting 
from an input for the parton distributions with A ~ 0.25 at values of Q 2 ~ 4GeV 2 was taking 
the best of both worlds. However, this lacked a real justification for the choice of input, which 
was significantly steeper than that expected from non-perturbative physics, but rather smaller 
than that from the BFKL equation, and also ignored potentially important ln(l/x) terms in the 
evolution. 

A significant step forward in the investigation of small-x structure functions was the develop- 
ment of the fey-factorization theorem [13] [14], the prescription for the way in which an off-shell 
photon-gluon scattering amplitude can be convoluted with the unintegrated gluon Green's function 
calculated using the BFKL equation to provide the small-x structure functions themselves. Hence, 
it enables one to find effective moment-space coefficient functions within the BFKL framework. 
Numerical calculations performed using this method were able to match the available data in a 
qualitative manner [15], as did similar calculations [16] using a modification of the BFKL equation, 
i.e. the CCFM equation [17]. The /c^-factorization formula was also shown to be very important 
if one insisted on working within the rigorous framework of the traditional renormalization group 
approach [18]. By showing how ^-factorization fits within the collinear factorization framework, 
Catani and Hautmann were able to calculate all the renormalization group anomalous dimensions 
to lowest nontrivial order in a s for each power of ln(l/x), and similarly for a number of coefficient 
functions. They demonstrated that within the renormalization group framework the Mellin trans- 
formations of the splitting functions, the anomalous dimensions "f gg (N,Q 2 ) and 7^ (TV ,Q 2 ) were 
the same renormalization-scheme-independent expressions as the effective anomalous dimensions 
[19] given by the BFKL equation, i.e. 

7 5 ° 9 (a s (Q 2 )//V) = f; aoJ^p-Y, ^ gf (a s (Q 2 )/N) = ^° gg (a s (Q 2 ) /N), (1.1) 

771=1 ^ ' 

where ^y gg (a s (Q 2 )/N) is given by the iterative solution of 

1 = ^Px(l° gg ), X(7) = 2V(1) " V>(7) " ^(1 " 7)- (1.2) 

and the solution as a power series in a s (Q 2 )/N exists for \N\ > \(a s (Q 2 )). They also derived 
expressions for ^j{a s (Q 2 )/N) and jj g (a s (Q 2 )/N) (series of the same form as in (1.1) but with 
an extra power of a s ) in certain factorization schemes (7° f{a s (Q 2 ) /TV) and ^ g (a s {Q 2 )/N) being 
zero), and also for the coefficient functions C 9 L1 {a s {Q 2 ) /TV), C f L1 {a s {Q 2 )/N), Cf^a^Q 2 )/^) 
and C{ 1 (a s (Q 2 )/N) (all zeroth-order quantities being zero except C^q, which is unity). This 
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facilitated calculations of structure functions within the normal renormalization group framework, 
but including much of what is often called BFKL physics (i.e. the leading-ln(l/x) terms), and 
indeed, a number of calculations were performed using somewhat different approaches [20]-[23], 
and in most cases comparisons with data made. It seemed that by including these terms it was 
not possible to improve upon the best fits for the small-x data using one- or two-loop evolution 
from soft inputs [20] [22]. Indeed, many ways of including them made the fits significantly worse, 
and this seemed to be universally true if the fits were more global, i.e. constrained by large-x data 
[20]. Also, it seemed that there was a very strong dependence on the factorization scheme used to 
perform the calculations when including the leading-ln(l/x) terms [20] [22] [24] [25], and a number 
of new factorization schemes were invented. 

The high precision of the most recent HERA data constrains theory far more than previously, 
and has changed the above picture somewhat. The best recent global fits seem to come from those 
intermediate approaches which use NLO perturbation theory with a quite steep (and completely 
unexplained) input for the singlet quark with A ~ 0.2 and a very Qf-sensitive small-x input for 
the gluon [26]. Fixed order perturbation theory using flat or valence-like inputs and low Q 2 j fails 
at the lowest x values, and relatively steep inputs for the singlet quark, i.e. A <S 0.2, seem to be 
absolutely necessary [27]. Approaches including the leading-ln(l/x) terms seemed to fail [28] [29] 
in practically all factorization schemes. 

In this paper we will take issue with the above conclusions. In particular we will demon- 
strate that the apparent failure of approaches using the leading-ln(l/x) terms, and certainly the 
factorization scheme dependence, is due to incorrect methods of incorporating these terms. We 
will show that the correct leading-order, renormalization-scheme-consistent (RSC) calculation of 
the structure functions is quite naturally factorization-scheme-independent, and also naturally in- 
cludes leading— In (1/x) terms in a well-defined form as well as terms leading-order in a s . We will 
also show that this calculation is clearly preferred by data. 2 

This paper will be structured as follows. We begin by defining conventions and giving a brief 
outline of the different possible types of scheme dependence in the calculation of structure functions. 
We then quickly review the work of Catani, who has already shown how to obtain factorization- 
scheme-invariant results in the small-x expansion by writing evolution equations in terms of the 
physical quantities, the structure functions and "physical anomalous dimensions", rather than in 
terms of parton densities and of the usual anomalous dimensions [31]. We then give a description of 
the calculation, within moment space, of structure functions in various expansion schemes using the 
normal parton language. The first part of this, regarding the normal loop expansion, will highlight 
points usually not discussed, particularly the form of the inputs. The second part, discussing 

2 We note that a very brief presentation of the complete RSC calculation of structure functions has 
already appeared in [30]. 
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the leading-ln(l/x) expansion, will present the only correct way to perform this expansion. One 
finds that by calculating physical quantities in a well-ordered manner the physical anomalous 
dimensions automatically appear, but not in the form one would obtain by simply solving the 
evolution equations using these anomalous dimensions at some given order. Moreover, there is 
a certain degree of predictive power for the form of the inputs for the structure functions at 
small x. We also explain why the standard solutions using the small-x expansions are strongly 
factorization scheme dependent. To conclude this section we present the argument that there is a 
unique renormalization-scheme-consistent calculation of structure functions, which applies to both 
large and small x. We present this calculation for both the currently academic case of the nonsinglet 
structure functions and for the phenomenologically important case of the singlet structure functions. 
We then briefly discuss how we move from moment space and obtain our x-space solutions, and 
the qualitative form these solutions must take, i.e. our best attempt at predictions. After this 
long theoretical presentation we consider the comparison with experiment, fitting the data for 
Tiix-, Q 2 ) using the renormalization-scheme-consistent solutions, and comparing to global fits using 
the normal loop expansion at NLO. We conclude that the full renormalization-scheme-consistent 
calculation gives the best global fit to structure function data, particularly at small x. Finally we 
investigate the phenomenological consequences for J : l(x,Q 2 ), and also preliminary indications for 
the charm structure function. 3 

2. Scheme and Scale Choices. 

For simplicity we work in moment-space for much of this paper, i.e. define the moment-space 
structure functions by the Mellin transform, i.e. 

F(N,Q 2 )= [ x N - 1 T{x,Q 2 )dx. (2.1) 
Jo 

The moment-space coefficient function is defined similarly but, as with the definition of the anoma- 
lous dimension, we define the moment space expression for the parton distribution as the Mellin 
transform of a rescaled parton density i.e 

f(N,Q 2 )= [ x N f{x,Q 2 )dx. (2.2) 
Jo 

The most general moment-space expression for a structure function is the sum of the products 
of the expressions for hard scattering with a certain parton (the coefficient functions) with the 
corresponding, intrinsically nonperturbative parton distributions, 4 

F(N, Q 2 ) = £ C a (N, a s G4), Q 2 / M |, /£/ »%)f a (N, a s (/x|), ^//4)> ( 2 -3) 

a 

3 There is a rather longer version of this paper which includes much more detail of a pedagogical nature 
and more comparison with other approaches [32] 

4 We deal with the heavy quark thresholds in a rather naive manner, i.e. the quarks are taken to be 
massless, with a particular flavour becoming active only above a certain Q 2 . 



4 



where \ip is a factorization scale separating the ultraviolet physics from the soft infrared physics. 
The coupling a s {^ 2 R ) satisfies the renormalization group equation, 



da(fi 



d\n(n 2 R ) 



£6„a"+V)^-/?(a(^)), (2.4) 



n=0 



where /j,r is the renormalization scale, the scale at which the coupling is defined. 

Two of the choices left open in the above expressions are common to all perturbative calcula- 
tions in quantum field theory: the choice of renormalization scheme (the choice of which expansion 
parameter is to be used, and consequently the form of the perturbative expansion), and subsequently 
the choice of renormalization scale. Conventionally the MS renormalization scheme is used, along 
with the simple choice [Ir = ^f, and indeed, the data seem to favour this choice of scale [33]. 
Hence, we take this simple choice. The remaining ambiguities are due to the particular problems in 
calculating quantities in QCD, i.e separating the physical quantity into the perturbative coefficient 
function and the intrinsically nonperturbative parton distribution. We have the freedom of choos- 
ing the factorization scale fip- As with renormalization scheme dependence this does not affect the 
all-orders calculation and invariance under this choice of scale leads to the evolution equation 

^| 1 = E^(« S (^))/ b (^). (2-5) 

It is desirable to choose \j? f to be large in order to make the expansion parameter a s (/j, F ) small, and 
to sum large logarithms in Q 2 , and fj, F is nearly always chosen to be equal to the hard scattering 
scale Q 2 . We shall also make this simple choice. 
This leaves us with our defining equations 

F(N,Q 2 ) = ^2c a (N,a s (Q 2 ))f a (N,Q 2 ) and = £ lab {a s (Q 2 ))f b {Q 2 ), (2.6) 

a b 

or being more careful there are two independent structure functions F^N, Q 2 ) and Fl(N, Q 2 )- In 
general we may write 



N 



Fi(N, Q 2 ) = -LfsTeAFfi^Q 2 ) + F? S (N,Q 2 ) (i = 2,L), (2.7) 

/ \ j=1 J 

where the singlet and nonsinglet structure functions are defined by 

if (AT, Q 2 ) = C{(N, a s )f s (N, Q 2 ) + Cf(N, a a )g(N, Q 2 ), (2.8) 

N f 

F? S (N, Q 2 ) = C? S (N, a.) £ ejf™{N, Q% (2.9) 



where Nf is the number of active quark flavours, f s (N,Q 2 ) and f^ s {N,Q 2 ) are the singlet and 
nonsinglet quark distribution functions respectively, and g(N,Q 2 ) is the gluon distribution. The 
equations for the nonsinglet distributions are ordinary differential equations, 

d f NS (N O 2 ) 

\l Q 2 =^s(N,a s )f- s (N,Q 2 ), (2.10) 
while those for the singlet sector are coupled 



d (f S (N,Q 2 )\ = fj ff (N,a s ) lfg (N,a s )\ ( f S (N,Q 2 ) 
dlnQ 2 V 9(N,Q 2 ) J \ lgf (N,a s ) j gg (N, a s ) ) \ g(N, Q 2 ) 



(2.11) 



This still leaves us one more freedom in our calculation, i.e. how we choose to remove the 
infrared divergences from the bare coefficient functions and hence how we define our parton dis- 
tributions. Starting from any particular choice for the definition of parton distributions one may 
always choose a new set of parton distributions by an invertible transformation 

f a (N, Q 2 ) = U ab (N, a s )f b (N, Q 2 ), (2.12) 

b 

where U a b(N, a s ) has a power series expansion in a s such that U a b{N, a s ) = 5 a b + 0{a s ). The struc- 
ture functions will clearly be unchanged as long as the coefficient functions obey the transformation 
rule 

C a (N,a s ) = (t/ T )- 1 (iV,« s )C 6 (iV,a s ). (2.13) 

By substituting f a (N,Q 2 ) = ^2 b U~ b 1 (N,a s )f b (N,Q 2 ) into (2.5) we easily find that the new par- 
ton densities evolve according to the standard evolution equations but with the new anomalous 
dimensions 

% b (N, a s ) = J2J2 U *^ N > a shcd(N, a s )U^(N, a.) + ^ (3{a s ) " s) U^N, a,). (2.14) 

c d c s 

The transformation defined above is called a change of factorization scheme. 5 However, it is 
important to realize that, unlike the changes in renormalization scheme, a change in factorization 
scheme leaves the expression for the structure functions unchanged not only to all orders, but order 
by order in a s . Calculations performed carefully at a given well-defined order in one scheme will 
lead to precisely the same results for the structure functions as those in another scheme. Also, the 
coupling constant to be used depends only on the ultraviolet renormalization. We will illustrate 
this point in §4. However, we will first discuss Catani's recent proposal for the construction of 
factorization-scheme-independent structure functions. 



5 The matrix U must obey a number of conditions in order that physical requirements on the parton 
distributions are maintained, e.g. flavour and charge conjugation invariance, fermion number conservation 
and longitudinal momentum conservation (see for example the second of [18]). However, none of these 
needs to be satisfied simply in order to keep the structure functions unchanged. 
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3. Evolution Equations for Structure Functions. 

It is, as Catani noticed [31], very simple to obtain factorization-scheme-independent 
factorization-scheme-independent effective anomalous dimensions governing the evolution of the 
structure functions In order to obtain these effective anomalous dimensions all one has to do is 
eliminate the parton densities from the equations (2.6) and (2.5). In order to demonstrate this, 
let us first consider the simple case of the nonsinglet structure function F 2 NS (N,Q 2 ). Multiplying 
both sides of (2.10) by YlJ=i e ] we can clearly write 



d 



{cM N . a .)) = ™s(N, ail ) (3.1) 



dh\Q 2 V C? S (N, a s )J '»*^>"°> C? S (N, a 
which becomes the factorization-scheme-independent equation 

dFj">{N,Q*) n _ . „ s 



dlnQ 2 



T 2 , NS (N,a s )F 2 Nb (N,Q 2 ), (3.2) 



where T 2yN s(N, a s ) = jNsiN, a s ) + d \n(C 2 NS (N, a s )) / d\n.Q 2 . Therefore, we have an effective 
anomalous dimension governing the evolution of each of the nonsinglet structure functions, and 
clearly T^ S (N, a s ) must be a factorization-scheme-independent quantity (and is in principle mea- 
surable). The solution to this equation is trivial: 



F 2 NS (N,Q 2 ) = F 2 NS (N,Q 2 )e W 



f-\nQ 2 

/ T 2tNS (N,a s )d\nq 2 



(3.3) 



The situation for the singlet structure functions is more complicated. However, using (2.8) 
for i = 2, L we may solve for the parton densities in terms of the two structure functions and the 
coefficient functions. Substituting these into (2.11) we then obtain the coupled evolution equations 

d (Fj(N,Q 2 )\ _ (T 22 (N,a s ) f 2L (N,a s ) \ (F 2 S (N,Q 2 )\ 
dlnQ 2 \F£(N,Q 2 ))- Vf L2 (iV, as ) t LL (N,a s )J \F*(N,Q 2 ) ) " [6A) 

The expressions for the physical anomalous dimensions, T22{N,a s ), T2L(N,a s ), etc., are straight- 
forward to derive in terms of anomalous dimensions and coefficient functions in any particular 
factorization scheme using the above procedure, but result in rather cumbersome expressions. It is 
simplest first to define a factorization scheme such that F 2 (N, Q 2 ) = f s {N, Q 2 ), i.e. C 2 (N, a s ) = 1, 
C 2 (N,a s ) = (this is generally known as a DIS type scheme [34] 6 ). In terms of the coefficient 



We call it a DIS "type" scheme because satisfying the above requirement still leaves freedom in how 
we may define the gluon density, and thus we are still considering a family of schemes. 



functions and anomalous dimensions in this type of scheme we have 

[(N,a s 



f 22 (N,a s ) = 7// (AT, a s ) - ^|^| 7/g (iV, a s ) , 



1 2L(JV,a s J — 



C g L {N,a s y 



T L , 2 (N,a s ) = C 9 L (N,a s ) lgf (N,a s ) - c£(N,a s ) lgg (N,a s ) + + C f L (N,a a ) Vf (N, aa ) 

- Ci(iV ' as) dln^ C£(iv>X 7/9(iV ' as) ' 

" / t» t- \ ,„ v dln(C?(AT,a s )) C((N,a a ) /Ar , 

W»J = ^< W '°'» + ding + £t^ 7/a(iV ' Q * ) - 

(3.5) 

Before going any further let us remark that we believe there is a (purely technical) problem with 
the above expression. As is well known Fl(N, Q 2 ) starts at an order of a s higher than F 2 (N, Q 2 ). 
Because of this there is an intrinsic asymmetry in the above definitions, with f 2 l(N, a s ) beginning 
at zeroth order in a s , r 2 2(N, a s ) and Tll(N, a s ) beginning at first order, and Tl 2 (N, o s ) beginning 
at second order. This asymmetry leads to the results obtained by solving these equations order 
by order being different that those using the normal parton picture even in the loop expansion. A 
trivial modification of Catani's approach is therefore to accept that Fl(N,Q 2 ) contains an extra 
power of q s , and to define the new structure function Fl(N,Q 2 ) = Fl(N, Q 2 )/(a s /(27r)). The 
longitudinal coefficient functions are likewise changed to C^(N,a s ) = Cf i (N,a s )/(a B /(27r)), and 
the singlet evolution equations become 

_±_ (F 2 S (N,Q 2 )\ _ (T 22 (N,a s ) T 2L (N,a s )\ (F 2 S (N,Q 2 )\ 
dlnQ 2 \F£(N,Q 2 )) ~ \r L2 (N,a s ) T LL (N,a s ) J {fE(N,Q 2 ) J ' ^ 

where the T(N, a s )'s are defined precisely as in (3.5), but in terms of of C1(N,a s ) rather than 
C1(N,a s ). This procedure restores the symmetry between the physical anomalous dimensions, 
and makes the order-by-order-in-a s calculations essentially the same as when using evolution of 
parton distributions. It is, of course, trivial to obtain the physical Fl(N,Q 2 ) from Fl(N,Q 2 ). 7 

Having made our redefinition of the quantities with which we work, we now have a direct 
relationship between possible calculations using the evolution equations for structure functions and 
the solutions using the parton densities. At present the parton anomalous dimensions and coefficient 
functions are known to order a 2 . It is easy to see that this allows us to derive each of the T's to 
order a 2 . Similarly, from the known expansions of the parton anomalous dimensions and coefficient 
functions in the form a™ X}m=i-n a m (a s /N) m [18], we can calculate T^^N, a s ) and r^ 2 (A^, a s ), 
T2 L {N,a s ) and F^N, a s ) (where both are trivially zero), and Tl L (N,a s ) and T^iN, a s ). This is 



Similarly, it is also best to work with the rescaled nonsinglet structure function F L (AT, Q ) . 
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the same order as for the parton anomalous dimensions, with the longitudinal anomalous dimensions 
having similar structure to the gluon anomalous dimensions and the r 2a (A, a s )'s having similar 
structure to the quark anomalous dimensions. The resulting expressions are relatively simple, being 

1 oo 

TlMs/N) = -j^rdcUv - tCl 1>0 )(f 7 ° 5 (a s MO + £(7°>*M0) n ) 
-frfi{N,a.), 

rl L (a s /N) = -J-a^ajN) + f>°> s /iV)n, ^ 

^ * n=0 

Tl 2 (N,a s ) = -{d f Lili0 - ICl'D^ajN), 

Tl L (a s /N) = 1 ° gg (a s /N), 

where 7V , C[ 10 and C 9 L 1 are the one-loop contributions to 'yj g (a s /N), C[ 1 (a s /N) and 
C 9 L ^atg/N) respectively. Each of these anomalous dimensions is renormalization scheme invariant 
as well as factorization scheme invariant, as we would expect for leading-order physical quantities. 

Solving the evolution equations for the structure functions using any subset of the currently 
known physical anomalous dimensions guarantees a result which is factorization scheme indepen- 
dent. However, there is considerable freedom in how we may solve the equations. We could, for 
example, simply put all of the anomalous dimensions currently known into (3.6) and then find the 
whole solution. Alternatively, we could solve using just the order a s anomalous dimensions and 
then try to perturb about this solution in an ordered manner. These two approaches would lead 
to rather different answers, but both would be factorization scheme independent. The problem of 
obtaining a correctly ordered solution for the structure function will be discussed in detail in the 
next section. We will initially use the familiar parton distributions and coefficient functions, and 
show that, even when using this approach, if we solve producing a well-defined expansion for the 
structure functions, we constrain the form of the solution rigidly and thus automatically avoid the 
problem of factorization scheme dependence. In fact, not only do we obtain the physical anomalous 
dimensions, but also find precisely how we should use them. 



4. Ordered Calculations of Structure Functions. 

There are in principle many different expansion methods one may use when obtaining solutions 
for the structure functions. The standard one is simply solving order by order in a s . But there is 
also the expansion in leading powers of ln(l/x) for given powers in a s (or equivalently in powers 
of N~ x in moment space), as we have already discussed. One can also combine the two expansion 
methods, and indeed, we will later argue that this is the only correct thing to do. Nevertheless, we 
will begin by outlining the procedure for making a well-ordered calculation of structure functions 
using the standard loop expansion. Although this is well known, we feel it is worth presenting it 
briefly, and making some points which are not usually highlighted, especially concerning the role 
of the starting scale. We may then discuss the more complicated cases of the expansion in leading 
powers of ln(l/x) and the combined expansion. 
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4.1. Loop Expansion. 

We begin by introducing some new notation. In order to solve the evolution equations for the 
parton densities order by order in a s and hence obtain expressions for F i NS (N, Q 2 ) and Ff(N, Q 2 ) 
we make use of equation (2.4) to rewrite the evolution equation for the nonsinglet parton density 
as 



» 2 s(Q 2 ) L , , ' 2 ? = -jNs(N,a s )f» s (N,Q 2 ), (4.1) 



d f£ s (N, Q* 
da s (Q 2 

where 7Ars(iV, a s ) = a 2 (Q )jns(N, a s )/ (3(a s ), and similarly for the singlet distributions, with 
similar definitions for jff(N, a s ) etc. as for 7jvs(-^V, ct s ). Each of the 7's may now be written as 



^N,a s ) = Y,r>\N)a: +1 , (4.2) 

n=0 

where where the super-subscript n, I denotes the (n + l)-loop quantity and we have an analogous 
definition for the normal anomalous dimensions. 

Using our definition, the evolution equations may be solved order by order in a s . When 
doing this it is necessary to choose a starting scale Q 2 for the perturbative evolution of the parton 
distributions (or equivalently a starting value of the coupling a s (Q 2 )), and specify input parton 
distributions at this scale. Let us discuss the choice of this scale briefly. Q 2 must clearly be chosen 
to be large enough that perturbative evolution should be reliable, i.e. a s (Q 2 ) 5s 0.3 and also such 
that "higher twist" (-^qcd/Q 2 ) corrections should be very small. Traditionally evolution only 
takes place up from this starting scale, for the simple reason of convenience, and also because some 
form of the inputs has been expected at low starting scales. Until a few years ago the requirements 
described above led to a choice of Q 2 ~ 4GeV 2 . In the past couple of years this value has often been 
chosen to be rather lower, due to the apparent success of evolution from lower starting scales, and 
also because much of the interesting small-x data is now at Q 2 < 4GeV 2 . These choices have been 
accompanied by guesses for the form of the inputs at low starting scales, e.g valence-like [9], or flat 
[5] [10]. We make no assumptions of the above sort about the value of Q 2 . We require it to be high 
enough to be in the perturbative regime and to avoid higher twists, but acknowledge that there is 
no reason why Q 2 cannot be chosen to be quite large, and evolution away from the starting scale 
performed both up and down in Q 2 . Taking this open-minded approach we also assume that only 
perturbative effects can lead to deviations from soft behaviour of the structure functions, and also 
demand that the form of our well-ordered expressions for the structure functions is as insensitive 
as possible to this choice, thus making the choice of Q 2 as open as possible. 

We begin by solving for the nonsinglet parton distributions, which are an easily understandable 
model. In this case the solution is particularly simple. Integrating both sides of (4.1) we obtain 



f£ s (N,Q 2 ) 



k=0 



2^ a s(Qi)J qj ,k( N iQi)\ \ as (Q2) ) exp 2^ ^ns( n > 



n 

L n = l 



(4.3) 
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Perhaps unconventionally, we explicitly express the input f^. s (N,Qj) as a power series in a s (Qj). 
This is necessary because changes in the starting scale Qj lead to a s (Qf )-dependent changes in the 
expression for the evolution term, which must be compensated for by a s (Qj)-dependent changes 
in the starting distribution in order to leave the whole expression for the parton distributions 
unchanged, as required. Let us examine this briefly by looking at the change of the lowest-order 
piece of (4.3), i.e. 

&,Q 2 ) = Co(iV,Qf)(^i|) 7Ws(JV) , (4.4) 

under a change in starting scale, Qj — > (1 + S)Qj, where 5 is some constant. We may regain 
expressions in terms of Qj by expanding the coupling constant in the form 

«s((l + 6)Qj) = a s (Qj) - 5b a 2 s (Qj) + 0(a 3 s (Qj)). (4.5) 
Under this change in the input coupling constant, the evolution term in (4.4) undergoes a change 

A wm) = - as{Qi)6b °^ s{N \^^) +higher order in asm - (4 - 6) 

This change in the parton distribution due to the variation in the leading term can be countered, 
up to higher orders, by a change in the order a s (Qj) input of the form 

A/£,!(AT, Qj) = 6b ^ l s (N)f» s (N, Qj). (4.7) 

Hence, changes in the evolution term due to a change in Qj begin at first order in a s (Qj), and 
are therefore absorbed by terms in the input at first order and beyond. The zeroth-order input is 
insensitive to such changes and is Qf-independent: f^^(N,Qj) = f^^(N)- Higher-order changes 
in the parton distributions due to changes in Qj are all accounted for by changes in the higher- 
order inputs which are equal to functions of N dependent only on the anomalous dimensions. 
Therefore, these higher-order inputs consist of perturbative parts multiplying the fundamentally 
nonperturbative f^fl(N). As we will soon discuss, there are other constraints to be satisfied, 
e.g. factorization scheme independence of the input for the structure function, and this slightly 
complicates the above picture, but does not change the main conclusions. 

Thus, it is necessary to express the input as a power series in a s (Qj) in order to make the 
parton distribution Q|-independent 8 , but only one intrinsically nonperturbative input which is 

8 Another reason for explicitly writing the input as a power series in a s (Qj) is that it makes little 
sense to demand that the starting distribution should not have a perturbative expansion, unless there is 
something special about a particular factorization scheme. Any transformation from a scheme in which 
the starting parton distribution is purely zeroth order in a s (Qj) will lead to a nonperturbative starting 
distribution multiplied by a power series expansion in a s (Qj). 
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(^-independent is needed. Usually in analyses of structure functions the parton inputs are taken 
to be a single a s (Q 2 )-independent function which is implicitly allowed to be Qf- dependent. Phe- 
nomenologically, this is normally much the same, but we stress the formally correct expression for 
the input here since it is rather important when constructing properly ordered solutions, and leads 
to some predictive power, especially in the small x limit. 

Substituting our solution for the parton distribution into (2.9), we obtain the general expression 
for the nonsinglet structure functions, 



F»»(N,Q 2 ) 



N 



J oo 



+ E Ci^xw 2 )) E e ? E ak M)& N > Qf) 



m=l j = l k = 



a s {Q 2 



cxp 



^ «(Q|)-<(Q 2 )) ^ - 



n 

L ri = l 



(4.8) 



It is clear that this may be written as 

/ (n2\\7Ns( N ) 00 n 00 

k ns (n,q 2 ) = E E ^MQ^r-iQ^TiQ 2 ) = J2 f ^ s (n,q 2 ), 

\<*sW ) J „ =0 m=0 n=0 

(4.9) 

and that, once a choice of renormalization scheme and starting scale have been made, each of the 
Fnm(N) must be invariant quantities under changes of factorization scheme. Any well-ordered 
calculation of the structure function should include all complete terms in (4.9) up to a given order 
in n and m, and no partial terms. In practice it is possible to work to a given order in n including 
all m < n, i.e. to expand to a given order in powers of a s (Q 2 ) plus powers of a s (Q 2 I ), if the 7's 
and Cj's are known to this order. 

We shall briefly describe how to construct this ordered solution for the structure functions, 
working up from zeroth order. Consider first calculating F^^iV, Q 2 ) by working to zeroth order in 
C^ S (N, a s ), Jns(N, a s ) and the starting distribution (remembering that this is Qf- independent). 
To this order 

Using the one-loop expression for the running coupling, as is appropriate for a lowest-order cal- 
culation, each of the quantities in this expression is factorization scheme independent and also 
renormalization scheme independent. Therefore we have a consistent leading-order (LO) expres- 
sion. If we calculate F^ s (N,Q 2 ) the zeroth-order coefficient function is zero. However, the only 
contribution for n = 1 in (4.9) comes from working to first order in C^ S (N, a s ) and to zeroth order 
in jNs(N,a s ) and the starting distribution. This leads to the LO expression for F^ S (N,Q 2 ) of 

^f(^,Q 2 ) = « s (Q 2 )^f^iV)^ e |/^ (iV)^^J . (4.11) 
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Again, using the one-loop expression for the running coupling, every term in this expression is both 
factorization scheme and renormalization scheme independent, giving a well-defined LO expression. 

We now consider the first correction to these expressions. The first-order expression for the 
renormalization group equation is 



d f^jN^ 2 ) 



a ^ = -7^ s (iV,«, S )/^(iV,Q 2 ) - a s (Q 2 )7^ 5 (iV)/^o(iV,Q 2 ), (4.12) 

with solution 



a s (Q 2 )) 



&N,Q 2 )= (a s (Q?)-a s (Q 2 ))7^(^ 

J y ^ " (4.13) 

Multiplying /^o by a s (Q 2 )C^ v 1 s i (A r ) and adding to /^i, we clearly obtain all terms in the expres- 
sion (4.9) for F2 (N,Q 2 )) at n = 1. Adding this to (4.10) we obtain the factorization-scheme- 
independent expression for FJ? S (N, Q 2 ) up to n = 1. We note that this is not the same as finding 
the complete solution to the renormalization group equation including all terms in the anomalous 
dimension up to first order in 7 and multiplying the solution by the coefficient function up to first 
order. This procedure would involve the exponentiation of the anomalous dimension, and thus 
would include incomplete parts of the F^(N) 's for n > 1, and would be a factorization-scheme- 
dependent, and hence physically ambiguous quantity. 

Similarly to F^ s (N,Q 2 ) we can obtain the NLO factorization-scheme-independent expres- 
sion for F^ S (N, Q 2 ). The expression at n = 2 is obtained by adding the product of the first-order 
coefficient function and f^ S i(N,Q 2 ) to the product of the second-order coefficient function and 
f£ S {N,Q 2 ). The NLO F^ S (N,Q 2 ) is the sum of this and (4.11). When working to NLO for 
either structure function we now have expressions which are renormalization scheme dependent. 
This scheme dependence compensates for the renormalization scheme dependence of the two-loop 
coupling constant (which has a renormalization-scheme-dependent value for Aqcd), and it is this 
expression for the coupling that we should use at this level. Doing so guarantees the renormaliza- 
tion scheme independence of the structure functions up to corrections of higher order in a s , i.e. 

C( a s-^2(L),0(l))- 

It is now simple to see how to construct structure functions order by order. Defining the n t h 
order renormalization group equation by 



d f g %(N,Q 2 ^ 
d*s(Q 2 ) 



a s (Q 2 ) JQ Z = " E l7s(N,a, S )f£ s n _ m (N,Q 2 ), (4.14) 



it is easy to prove by induction that the solution contains all terms in the full solution with a given 
sum of powers of a s (Q 2 ) and a s (Qj) multiplying the everpresent (a s (Qj)/a s (Q 2 )) ;:f Ns( N ) factor. 
Thus, F" n s (N,Q 2 ) is given by 

Nf n 

Fl? n S (N, Q 2 ) = £ e) £ C?^a™{Q 2 )fl%_ m {N, Q 2 ). (4.15) 

j m=0 
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The nth-order scheme-independent structure function is then the sum of the Fj^f(N,Q 2 ) up to 
order n. Including all FJ?J^(N, Q 2 ) up to order n is working to (n+l) t h nontrivial order, and requires 
the (n + l)-loop coupling in order to make the expression renormalization scheme invariant up to 
higher orders in a s . Similarly, including all F^j^ (N, Q 2 ) up to order n is working to n t h nontrivial 
order, and requires the n-loop coupling. 

This procedure clearly provides factorization scheme independence and renormalization scheme 
independence for this method of expansion. We can also see how it relates to the discussion of the 
factorization-scheme-invariant evolution equations in terms of the structure functions. In order to 
do this let us consider the solution for the non-singlet structure function F^ S (N,Q 2 ) again. We 
may rewrite our general solution (4.8) in the form 



N f 

?N S ( at 



F?»(N,Q') = 



m=l j = l k=0 

r 00 („,n(n2\ „,n(r>2\\ />ln 



a s (Q 2 ^^ 



exp 



oc s {Q 2 

l„ =1 n J\nQ] amq m=i 



X 



(4.16) 

This way of writing F 2 NS (N, Q 2 ) is particularly useful since it separates the solution into the value 
of the structure function at Q 2 (the term in square brackets), and the ratio of its value at a 
different Q 2 with this initial value (the rest of the expression) . Clearly these two quantities must 
be separately factorization scheme independent. Also, it is clear that this solution is of exactly the 
same form as (3.3), and we can express it simply in terms of an input for F 2 NS (N,Q 2 ) at Q 2 and 
a physical anomalous dimension which governs the evolution, i.e. our input and our prediction. 

The input and the evolution will mix with each other if we make a change in the starting 
scale, however. Examining the effects of such a change for the full physical quantity gives us 
information about the form of the input. The lowest-order input for the structure function is 
just F^q(N) = Ylf=i e2 jfqj%{N). Making the change of starting scale and consequently of a s (Qj) 
already considered, the change in the lowest-order evolution is as in (4.6), and this leads to a 
change in the lowest-order structure function which is of higher order, and which can be absorbed 
by a change in the NLO input for the structure function of 

AF 2 N f(N, Q 2 ) = Sb ^ l s (N)F 2 N s (N). (4.17) 

In terms of parton distributions and coefficient functions 

F^fiN, Qj) = £ e 2 {f^(N, Q 2 ) + C 2 ^(iV)/^ (iV)). (4.18) 
.7=1 

We chose the change in f^(N,Qj) in (4.7) so that the structure function would be independent 
of Qj, and it is clear that that is consistent with (4.17) and (4.18). However, we can now say 
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more about the form of f^ S i{N, Qj). Because it is the leading term in the expression for the input 
which depends on a s (Qj), F 2 N f(N,Qj) must be renormalization scheme independent. However, 
C^fj(N) is renormalization scheme dependent, so f^{N, Qj) must also be renormalization scheme 
dependent in a way such as to cancel this. Hence, fq- S \{N,Qj) must not only have a part like 
ln(Qf )7^(AT)/^(AT) in order to maintain (^-independence of the structure function, but also a 
part like —C^f^N)/^^) in order to maintain renormalization scheme independence, i.e. 

&N, Qj) = (HQj/A NS )^ l s (N) - CSf tl (N))fSfim (4-19) 

where A^s is some unknown scale parameter. So we see that 

F 2 N f(N, Qj) = ln(Ql/A NS )^ l s (N)F 2 N s (N) ee ln(Qj/A NS )T% l s (N)F 2 N s (N). (4.20) 

It is clear that this does not spoil our argument that f^^(N, Qj) consists of perturbatively calcu- 
lable quantities multiplying the fundamentally nonperturbative f^fi{N), and that all the higher- 
order inputs may be chosen to be perturbative functions multiplying this nonperturbative input, 
and therefore that the input for the structure function is a perturbative power series (depend- 
ing on the physical anomalous dimension) multiplying the single nonperturbative factor f^^{N), 
which may also be interpreted as a fundamentally nonperturbative input for the structure function 
F 2 q(N). Hence, demanding invariance of our expression for the structure function under changes 
in Qj leads us to a power series expression for the input, but with only one (for each quark) really 
nonperturbative factor for this input. We also see that if Qj = Ajvs, the first-order perturbative 
correction to F 2 NS (N,Qj) vanishes. Hence, we might expect A^s to be some scale typical of the 
transition between perturbative and nonperturbative physics, i.e A^s & lGeV. 

Separating the expression for the structure function into a definite input and evolution part 
also enables us to view the loop expansion in an alternative manner. We see that when expanding 
out to n t h order in the loop expansion we are including all terms where the order of the input 
part for the structure function added to the order of the evolution part of the structure function 
is less than or equal to n. Writing the solution as in (4.16) does, however, also illustrate that 
demanding factorization scheme invariance does not on its own force us into the strictly defined 
loop expansion. It is clear that we could, if we wished, expand the input and evolution term out 
to different orders in a s , still maintaining factorization scheme independence. However, this is not 
a sensible approach for reasons of renormalization scheme dependence. If we were to expand out 
the input and evolution terms to different orders in a s we should really use a s itself calculated to 
a different order in each case, surely a perverse thing to do. It makes far more sense to combine 
physical quantities to a given order. 

The solution for the longitudinal structure function is much the same as for F 2 N q(N, Q 2 ). The 
situation for the singlet structure function is Exactly analogous to that for the nonsinglet, but 
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is complicated by the fact that we now have coupled evolution equations for the quark and gluon 
distributions. This makes it impossible to write a closed form for the solution to the renormalization 
group equations in the way we did for the nonsinglet case in (4.16), but the equations may be solved 
order by order in the same way, i.e. the lowest-order solution is 



f b (N,Q 2 ) = fS' + (N)l +fS'-(N) 



a s (Q 2 )J U V 'WQ 2 

9o(N, Q )=g (N) j + g (AO j 



(4.21) 



where j+ l (N) and 7°' (A 7 ") are the eigenvalues of the zeroth-order matrix for the 7's, and fo' + (N) + 
fo'~(N) = fo(N) and go(N) + g$ (N) = go(N). As with the nonsinglet quark distributions, the 
lowest-order inputs for the partons are Qf— independent. Also, if we expect any powerlike behaviour 
to come only from perturbative effects, then these Qf-independent inputs for the quark and gluon 
are analytic for N > 0. 

We may also discuss the relationship to the solutions using the evolution equations for the 
singlet structure functions. Unlike the nonsinglet case it is rather difficult to see how to express the 
solution in terms of factorized inputs and evolution parts; the evolution parts will not be simple 
exponentials, as in (4.16). However, the solution for the structure functions can be written as 
a series of terms each of which can be factored into a part dependent only on a s (Q 2 ) and one 
depending on both a s (Q 2 ) and a s (Q 2 ), the latter either vanishing at Q 2 = Q 2 , or being equal to 
{a s {Q 2 )/a s {Q 2 )y+'-^~ 1 ^ and thus equal to unity at Q 2 = Q 2 . For each of these terms the former 
part may be interpreted as an input and the latter part may be interpreted as the evolution. Within 
the loop expansion we then include all terms where the order of the input part plus the order of 
the evolution part sums to less than or equal to some integer n. 

By examining the form of the inputs under a change in starting scale, as with the nonsinglet 
structure functions, we find that the only fundamentally nonperturbative parts of the inputs are the 
zeroth-order parts, F^ Q (N) and Ff (N), with all other inputs being in principle expressed as per- 
turbative functions of the physical anomalous dimensions multiplying one of these nonperturbative 
components. They are given in terms of the parton inputs by 

F 2 S o(N) = fi(N), (4.22) 

and 

Ko( N ) = CUN)f s (N) + Cl^g^N). (4.23) 

We choose to think of the expressions for the structure functions as the real nonperturbative 
functions, since they are, of course, the physically relevant quantities. 

The extra complexity of the solution in the singlet case, when compared to the nonsinglet 
case, also means it is far from obvious how to express the physical anomalous dimensions in terms 
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of the parton anomalous dimensions and coefficient functions simply by comparing the forms of 
the solutions in the partonic language and purely in terms of structure functions. However, the 
expressions for the physical anomalous dimensions can be found at each order using (3.5). One 
could then solve for the structure functions in a manner different from the loop expansion, e.g. 
calculate the solution to the whole evolution equation for the structure functions using the physical 
anomalous dimensions up to a given order. This would be by definition factorization scheme 
invariant, but would only be renormalization scheme invariant up to the same order in the solution 
as the order of the anomalous dimensions. The rest of the solution would contain only a subset of 
the possible terms of a given form obtained from the full calculation. This subset will not necessarily 
be a good representation of the full set of terms obtained at this order, and would therefore have 
no real significance, and should be dropped. Hence, the evolution equations in terms do nothing 
to alter the strict ordering of the solution using this method of expansion. 

This whole discussion of the order-by-order-in-a s expansion scheme is perhaps a little aca- 
demic since the errors invoked in performing the calculations without paying strict heed to the 
formally correct procedures are rather small. For example, when using this standard method of 
expansion only two different factorization schemes are generally considered, the MS scheme and the 
DIS scheme; the latter related to the former by a change of parton distributions, so that the singlet 
quark distribution is equal to the singlet structure function. If a calculation is made in one scheme 
to a well-defined order, and then a transformation to the other scheme made correctly, precisely 
the same result will be obtained. In practice, small differences are noticed between calculations 
using different schemes within the loop expansion to NLO, but these come from well-understood 
sources. One common source is that the starting distributions in both schemes are described by a 
simple functional form, e.g. 

xf(x, Qj) = A(Q 2 j)x- x( - Q2 i\l - x)^ Q '\l + eiQ^x 1 / 2 + j(Qj)x), (4.24) 

rather than the formally correct expression of a power series in a s (Qj) with essentially perturba- 
tive coefficients convoluted with a Qf-independent nonperturbative function of x. If the starting 
distribution is of the form above in one scheme, then in the other scheme it will not be modelled 
precisely by a function of the same form. However, the error is in general small. Alternatively, if 
the calculations are not done in a well-ordered manner, then differences between the calculation 
done in the two schemes, or between this type of calculation and the correct NLO calculation, will 
be of NNLO. Again this usually results in only small differences. 9 Similarly, small differences would 
be obtained by working in two different renormalization schemes. 

9 For a comparison of calculations done at NLO using different methods see [35]. 
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Nevertheless, there are a couple of points we wish to make here, concerning the loop expansion. 
Firstly we note that if one considers the input to be consistent with the loop expansion, e.g. NLO- 
order evolution should be accompanied by a NLO input, the power of A(Qj) should not correspond 
to parton distributions much steeper than flat for the singlet quark or gluon at this order. This is 
because the first-order-in-a s (<5/) input should be accompanied by no more than a single power 
of ln(l/x). Restricting X(Qj) in this way is rather important for the fits to the low x data, and 
would mean that NLO fits to small-x data would be very poor. The only way to avoid this is to 
let X(Qj) be an artificial free parameter, in which case, if (4.24) describes the singlet quark density, 
it must be ~ 0.2 — 0.3 for practically any Q\. This value is totally unjustified, and the need for 
this steepness in the input for the quark is a clear sign of the limited usefulness of the NLO-in-a<. 
calculation at small x. 

There is another reason for being concerned about the validity of the loop expansion at small 
x. The reason for the relative smallness of the differences between inaccurately performed NLO 
calculations noted above, even at small x, is that the differences between these calculations do not 
contain terms which are any more leading in ln(l/x) than the NLO calculation itself. Hence they 
are genuinely an order of a s down on the NLO calculation, with no small-x enhancement. This is 
also true for calculations done in different renormalization schemes (or at different renormalization 
scales). However, the real NNLO contribution are higher order in a s but also contain terms at 
higher order in ln(l/x), and are potentially large at small x. Hence, the relative insensitivity of 
structure functions to changes in renormalization or factorization scheme for calculations which are 
not carefully ordered is no guarantee that genuine higher-order corrections will be small at small 
x when using the loop expansion. Indeed we would naively expect them to be large. 

In contrast to the insensitivity when using the loop expansion, when using the leading-ln(l/x) 
expansion very large differences between calculations done in a large number of different factoriza- 
tion schemes have been noted. This is an clear sign that the calculations are not being done in 
a well-ordered manner, and that the ambiguity introduced by lack of care in the calculations is 
greater in this method of expansion than the standard loop expansion. We will now demonstrate 
that this is indeed the case. 

4.2. The Leading — ln(l/x) Expansion. 

It should clearly be possible to define a well-ordered expansion in leading powers of ln(l/x), 
or equivalently, in leading powers of 1/N in moment space. We will now demonstrate that this is 
indeed the case. As in the loop expansion, we will first work in terms of the traditional parton 
distribution functions and coefficient functions, and see how this results in expressions containing 
the physical anomalous dimensions. Doing this enables us to see how large factorization scheme 
dependence can arise when calculating less carefully within this expansion scheme. 
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First we must set up our notation. In this expansion scheme we write, 

oo oo oo 

7(iV,Q 2 ) = ^<(Q 2 ) Y, T m <(Q 2 )N- m ^Ya:(Q 2 )r(<x s (Q 2 )/N). (4.25) 

n=0 m=l — n n=0 

In particular, in the MS renormalization and factorization scheme we may write 

oo 

% g (N,a s (Q 2 )) = Y J <{Q 2 )T gg MQ 2 )/N), (4.26) 

77 = 

where the series expansion for 7° (a s /iV) is known to all orders (all the coefficients being positive). 
This expression for j gg (a s /N) is renormalization scheme independent. 7* (a a /JV) is of course 
renormalization scheme dependent, since it must change to absorb part of the effect of the 0(a 2 ) 
change in the coupling on 7° {a s /N). The renormalization scheme independence is also true for 
7°j(a s /A^), which obeys (1.1), and also for ^®j(a s /N). There is also a renormalization-scheme- 
independent relationship between 7^(a s /iV) and jj g (a s /N), which tells us that 

^/N) = ^j g (a s /N)-^, (4.27) 

where the second term in the brackets is the one-loop contribution to jj g (a s /N). It is also known 
that neither the nonsinglet anomalous dimension nor the nonsinglet coefficient function have poles 
at N = 0. Hence, the nonsinglet sector makes very little contribution to the structure function at 
small x, and as such we will ignore it for the remainder of this section. 

A general change in factorization scheme may be expressed by writing an element of the 
transformation matrix U as 

00 00 00 

U ab (N,a s (Q 2 )) = Y,^(Q 2 ) K b m aT(Q 2 )N- m = J2<(Q 2 )K b MQ 2 )/N), (4.28) 

77 = 777=1 — 77 77 = 

with condition on the U™ b m such that U obeys U a b = 5 a b + 0(a s ). A change of factorization 
scheme with U® b 7^ can introduce scheme dependence into the 7° b (a s /A^)'s, as we see from (2.14). 
Insisting that ^ g {a s /N), j g ^(a s /N), j^(a s /N) and ^® g (a s /N) are unaltered by scheme changes 
leads to the requirement, 

U° fg (N,a s ) = 0, U° ff (N,a s ) = l, 

00 (4 29) 

U Q gg (a s /N) = 1 + J2 U Q gg n (a s (Q 2 )/N)\ U° gf (a s /N) = l(U° gg (a s /N) - 1). 

77=1 

This requirement also preserves the relationship (4.27) between jjj(a s /N) and jj g (a s /N). For 
simplicity, and due to factorization scheme invariance of physical quantities, we will only consider 
factorization scheme changes away from MS scheme of the type (4.29). The U^ b {a s /N) for n > 
will have no restrictions. 
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Restricting ourselves to these schemes we may write 

oo oo oo 

Ci(N,a s ) = l + ^2<(Q 2 ) £ C ln,m a 7(Q 2 ) N ~ m = 1 + £ a "(Q 2 ) C ln( a s/N), (4.30) 



n=l m=l — n 

and all other coefficient functions as 



n = l 



C?(N,a s ) = J2<(Q 2 ) £ C i,n, m <(Q 2 )N- m = J2^(Q 2 )C9 n (a a /N). (4.31) 



n=l 



m=l — n 



n=l 



All the Cf n (a s /N) are both renormalization-scheme- and factorization-scheme-dependent quan- 
tities. Indeed, all of the Cf n m are renormalization scheme and factorization scheme dependent, 
except for the C£ n l _ n , which come from the one-loop longitudinal coefficient functions which, as 
we saw in the previous subsection, are totally scheme independent. There are also two renormaliza- 
tion and factorization scheme (with our restrictions) independent relationships between coefficient 
functions: 

CtMs/N) = ^(c^iaJN) - Cf ili0 ), (4.32) 

where the second term in brackets is the one-loop contribution to Cf ^(ag/iV), which is itself renor- 
malization scheme and factorization scheme dependent (being equal to (Nf/Qit) in MS scheme); 
and 



(4.33) 



where the second terms in the brackets are the one-loop contributions to C£ 1 (a s /A^) and 
C 9 L 1 (o; s /iV), both of which are renormalization and factorization scheme independent. 

Working in an arbitrary factorization scheme (up to the above restrictions) and using the 
general expressions for the 7's and coefficient functions, we may find expressions for the structure 
functions. The first step towards this is solving the renormalization group equations for the parton 
distributions. The lowest-order part of the equation is, 



d 



fi(N,Q 2 ) 

9o(N,Q 2 ) 



da s (Q 2 

This may easily be solved to give 
fo S (N,Q 2 ) = f S (N,Qj), 
go(N, Q 2 ) = (g (N, Q 2 ) + § f$ (N, Qj)) exp 





§ 7 ° 5 (a s /A0 T° gg (a a /N) 



fo S (N,Q 2 ) 

g (N,Q 2 ) 



(4.34) 



a * {Q2 >HMs{q 2 )/N) 



U«,(Q 2 ) 



da s (q 2 ) 



-lfo(N,Q'j). 



(4.35) 

This is analogous to the lowest-order solution within the loop expansion and contains two 
factors, one of which must appear in all the higher terms in the expansion: instead of 
{a s {Q 2 ) /a s (Q 2 )y i +'- ( " N \ corresponding to the two eigenvalues of ^ 0,l (N) in the loop expansion, we 
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have exp 



r s(Q 2 r) l ag (» s (l 2 )/N) , , ^ 



and 1, corresponding to the two eigenvalues of 7°(a s /iV) in 



the leading-ln(l/x) expansion, i.e. ^® g (a s /N) and 0. 

The fact that the leading order part of 7 9g and 7 S / is more leading than the leading part of 7/ 9 
and 7//, in this expansion scheme makes obtaining a well-ordered solution rather more complicated 
than in the loop-expansion. (The same is true when using the physical anomalous dimensions since 
T 2 L and T 2 2 both equal 0.) Hence, before using our solutions for the parton densities to construct 
expressions for the structure functions we will need to solve higher-order renormalization group 
equations in order to determine the general form of the solutions. This is done iteratively using 
evolution equation ordered analogously to (4.14), and in order to get useful solutions we must go 
out to n = oo, keeping only the most leading terms each time we iterate. With some work and a 
straightforward inductive proof we obtain 



f(N,Q 2 ) =a s {Q 2 



j} g (a s (Q 2 )/N) 
l° gg (a s (Q 2 )/N) 



(go(N) + lf*(N)+g (N,Q 2 ))x 



exp 



^° gg (a s (q 2 )/N) ^ gg {a s { q 2 )/N) + ^} g {a s {q 2 )/N) 



a s (Q 2 ) 



a 2 s {q 2 ) 



+ 



a s (q 2 ) 



da s (q 2 



a s (Qj) 
a s (Q 2 ) 



4 2Nj, 
'96^ 



+ fi(N, Q 2 ) + higher order in a s and/or N, 



(4.36) 



and, 
g(N,Q 2 



<9o(N) 
exp 

- lfo S (N) 



lfi(N) + ~g (N,Qj))x 

^ l%{a s {q 2 )/N) jl g (a s (q 2 )/N) + lj} g (a s (q 2 )/N) ^ 

a 3 (Q2) a>f{q 2 ) a s {q 2 ) 

a s (Qj 



(4.37) 



4 2Nf 

' 9 e^b^ 



+ higher order in a s and/or N. 



MQ 2 ) 

We may also write the expression for the rate of change of the quark distribution 



2 ,df s (N,Q 2 ^ 



da s (Q 2 



=a s (Q 2 )f fg (a s (Q 2 )/N)(g (N) + |/ s (iV) + g (N, Q 2 ))x 

^(Q 2 ^ gg (a s ( q 2 )/N) ^ g (a s ( q 2 )/N) + ^) g (a s (q 2 )/N) 



exp 



/ 

■J a 



la s (Q 2 )f*(N) 



a 2 (q 2 ) + 
2N f fa s (Q 2 ) 



a s (q 2 



da s {q 2 ) 



67t6 \ a s {Q 2 



+ higher order in a s and/or N. 

(4.38) 

This last expression will be important since at leading order a 2 (Q 2 ){dF 2 {N , Q 2 )/da s {Q 2 )) is di- 
rectly related to dF 2 (N,Q 2 )/d InQ 2 which we will wish to study as well as F 2 (N,Q 2 ). 
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under 



The input go(N,Qj) is chosen so that the change in exp 
a change in Q\ can be compensated for by a change in go(N, Q"j) up to corrections of higher order. 
Hence, the gluon input may be written as 

(go(N) + lf s (N)+~ 9o (N,Qj)) 

= (9o(N) + |/ 5 (iV))(l + £ ffo, m (^^) m ) exp[ln(Qf/A ra ) 7fl ° g (a s (Qf)/iV)], 

(4.39) 

where A gg is an unknown scale. The series J2m=i 9o,m{oi s {Q 2 i) / N) m is at yet undetermined, but is 
potentially renormalization and factorization scheme dependent. It will only be determined when 
we come to construct the structure functions themselves. In the same manner we can determine 
the general form of our input, ff(N, Qj), obtaining 

fi(N, Qj) = (g (N) + lf s (N) + g o{ N,Qj))aM) 7 ^^^ + fl(N,Qj), (4.40) 
where 

oo / ( f\2 \ \ 771 

f 1 S (N,Ql) = N^fU^f l1 ) (9 (N) + lf s (N)) (4.41) 

n=0 ^ ' 

is potentially renormalization and factorization scheme dependent, but not yet determined. 

We also consider the form of the iV-dependence of our inputs. If we take the point of view 
that any steep behaviour in the parton distributions only comes about due to perturbative effects, 
then we assume that /q(N) and go{N) are both soft, i.e. either flat or even valence-like when the 
transform to x -space is performed (or at most going like a finite, small power of ln(l/x)). This 
requires that they both be analytic for N > 0. As in the loop expansion, we may think of go(N), 
and /q(N) as fundamentally soft, nonperturbative parts of the input. 10 The parts multiplying 
these are then really determined by perturbation theory. Since we are meant to be expanding our 
solution for the structure functions in powers of both a s and N, it might be argued that we should 
expand /q(N) and go(N) in powers of N. We feel this is not really appropriate since it is the 
perturbative part of the solution for which we are able to solve, and thus which we are able to 
order correctly, and in the expressions for the structure functions the whole of the nonperturbative 
inputs should multiply the well-ordered perturbative parts of the solution. 

It is now possible to examine the form of the solutions for the structure functions. We can con- 
struct the leading part of the full solutions by combining our solutions for the parton distributions 



10 Of course, we should be trying to find fundamentally nonperturbative inputs for the structure functions 
rather than the partons, as discussed for the loop expansion. However, in this expansion scheme fg (N) 
and go(N) are trivially related to F 2t o(N) and Fl,o(N) as we will see soon. 
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with the zeroth- and first-order coefficient functions. This gives 



F L (N,Q 2 ) = a s {Q 2 )C 9 L!l {a s {Q 2 )/N){g {N) + § f*{N) + g (N, Qj)) 



cxp 



^% g {a s {q 2 )/N) | % g (a s (q 2 )/N) + ^) g (a s (q 2 )/N) 



+ 



+ a s (Q 2 )(Cl h0 -ICl h0 )f s (N) 



a s (Qj) 
a s (Q 2 ) 



a s (q 2 

4 2Nf 



da s (q 2 



+ higher order in a s and/or N, 

(4.42) 



where (4.33) has been used. Also 



F 2 (N,Q 2 ) = as (Q 2 )^^^ + Cf> s (Q 2 )/iV))MA0 + f/o^W + ~g (N, Q 2 ))x 



exp 



/ 

■J a 



a ^° gg (Mq 2 )/N) l] a {a s {q 2 )/N) + ^} g (a s ( q 2 )/N) 



a 2 s (q 2 ) 



+ 



a s (q 2 ) 



da s (q 2 



(4.43) 



+ f S (N) 



« S (Q 2 ) 



4 2JVj 

" 9 eT^J 



+ ff(N, Qj) + higher order in a s and/or N, 



and 



2 2 ^F 2 (iV,Q 2 ) 



da s (Q 2 ) 



(a s (Q 2 )j} g (a s (Q 2 )/N)+ 



a s (g 2 )Cf jl (« s (Q 2 )/iV)7° g (a s (Q 2 )/iV))(5o(iV) + |/,f W + go(N, Q 2 )) 
' r^)^ gg (a s (q 2 )/N) ^ g (a s {q 2 )/N) + ^) g ( as (q 2 )/N) 

J a 



exp 



a 2 s (q 2 ) 
a s (Qj) 



a s (q 2 



da s {q 2 ) 



9 6vr6 ' 



a s (Q 2 ) 



4 till 



+ higher order in a s and/or N, 



(4.44) 



where, we drop the superscript S for the structure functions for the rest of this subsection. 

Each of these expressions will be factorization scheme independent since each represents the 
expression for a physical quantity up to corrections of a form different from the terms explicitly 
appearing and which we have deemed to be higher order in our expansion scheme. However, they 
are not in themselves the leading order expressions. In order to derive these we begin with the 
longitudinal structure function. Making the definition 



■ {Q H 9 Mq 2 )lN) + ^) g (a s (q 2 )/N) 
« S (Q 2 ) a s (q 2 ) 



d 



(4.45) 

where ^f(Q 2 ,Qj) is factorization-scheme independent (as can be checked using the rules (2.12)- 
(2.14)) and ^i{Q 2 ,Q 2 ) = 1, we may factorize (4.42) completely into input and evolution parts, 
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i.e. 



2 /V 

F L (N,Q 2 ) = a s (Q 2 )-^(g (N) + |/ s (iV))exp 



/ 2? 2^ (« 



— ("i. sl r ; 1 

'a s (Q 2 ) 



C£ ;1 (a s (Q?)/iV) / ffo(AT, Q 2 -) , 



C£,i,o A ft,(JV) + | f S (N) 

4 2W L _ 



<ft(Q 2 ,Q'j)-ln( aAQ2l] 



+ « S (Q 2 ) ( 18 27 l Nf ) fo(N) (^|||y) 9 + higher order in a s and/or N. 

(4.46) 

It is now possible to attach direct physical significance to each of the factorization-scheme- 
independent pieces appearing in this expression. 

We first consider the inputs. Going to (4.43) for the moment we see that fo(N) is the only 
term in F 2 (N,Q 2 ) which is zeroth-order in our expansion scheme. As such it is the zeroth order 
input for F 2 (N, Q 2 ), and we may write 

fi(N) = F 2t0 (N), (4.47) 

and this is one of our two fundamentally nonperturbative inputs. Also, from (4.46) we see that 
the total a s (Q 2 )-independent input for Fl(N,Q 2 ) (once we have divided out a single power of 
a s (Q 2 )/{2^)) is 

^(<7o(A0 + £f${N)) = F Lfi (N). (4.48) 

This is therefore equal to our other fundamentally nonperturbative physical input We also 
make a similar definition for the part of the input for F L (N,Q 2 ) which is of the form 
£~=i a n (Qj){a s (Qi)/N) n multiplying (g (N) + ±f$(N)), or more correctly multiplying (F Lj0 (N)- 
((36 - 8N f )/27)F 2 fi(N)), to complete our definition of the input in (4.46) and write 



(4.49) 



~^{9o{N) + lf s (N) + go(N, Q 2 )) ( C £ * J 

= F Lfi {N) - ( 36 ~ 2 * Nf y 2 ,o(N) + F L , (N,Q 2 ) 

This whole expression must be both factorization scheme and renormalization scheme independent, 
facts which reveal information about the form of the gluon input. C 9 L 1 {N^a s {Q 2 I )) is both renor- 
malization and factorization-scheme dependent, so the scheme dependence of go(N,Qj) must be 
precisely so as to cancel this out, Hence, there is no reason for F l,o{oi s {Q 2 i ) / N) to depend on the 
leading-order longitudinal gluon coefficient function at all, and indeed, a natural choice seems to 
be that (1 + YZ=i 9o, m ( a s{Qj)/N) m ) is chosen equal to {C 9 LX JC 9 L1 (N, a s (Q 2 j ))) (it is difficult 
to see what else it could be chosen equal to), and therefore 

F L ,o(N) - ( m ~J Nf ) F *A N ) + ~Fl,o(N, Qf) 

/ /36-8/V\ \ (450) 

= \F Li0 (N) - [ 2? f jF 2 AN)j exp[ln(Q?M LL ) 7 « g (« s (Q 2 )/iV)]. 
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Hence, we have a prediction for the input for the longitudinal structure function at small x in terms 
of the nonperturbative inputs and some scale All (where All = A gg from the previous subsection). 
As with Ans earlier, All is the scale at which the input is equal to the nonperturbative input 
alone, and hence we would expect it to be typical of the scale where perturbation theory starts to 
break down. Qj is a completely free parameter. We have constructed the solution to be formally 
insensitive to Qj at leading order, but there is clearly some residual Q 2 -- dependence. Hence, there 
will also be some optimum Qj to choose as the starting scale. 

We may now examine the terms governing the evolution, exp 



f a.(Qj) 7g B (q.(g 3 )/*0 . , /„ 2. 
JaJQ 2 ) aa s{q 



IS 



the scheme-independent factor governing the small-x growth with Q 2 , and we make the definition 

r«.(Q?) ^(a s (q 2 )/N) 

'a s (Q 2 ) 



K(Q 2 ,Qn 



f 

J a. 



-da s (q 2 



(4.51) 



al(q 2 ) 

Of course, the other evolution factor resulting from the eigenvalues of the zeroth-order anomalous 
dimension was simply unity and as such $q (Q 2 ,Qj) does exist, but is implicitly zero. However, 
there is a correction to this factor of unity in (4.42)-(4.44), and we make the definition 



*I(Q\Q 2 i) 



4 2N f lQ fa s (Q 2 



a s (Q 2 ) 



(4.52) 
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Having made these factorization-scheme-invariant definitions for the inputs and evolution we 
may write the solution for Fl(N,Q 2 ) as 

a s (Qj) 



F L (N,Q 2 



2tt 



-(F L , (N)+F Lfi (N,Qj) 



exp 



+ 



<s>t{Q 2 ,Q 2 i) + <s>t{Q 2 ,Q 2 i) 

36 - 8N f ' 



In 



a s (Qj 



2vr 



27 



F 2 ,o(iV)exp 



36 - 8Nf 
27 

* s (Qj) \ 
u s {Q 2 ))_ 

^(Q 2 ,Qj)- 



In 



a s (Q 2 ) 



(4.53) 



+ higher order in a s and/or N. 

This is still of mixed order, but it is now relatively obvious how we may separate out the "leading 
part" from this expression. [&f(Q 2 , Qj) — \n(a s (Qj) / 'a s (Q 2 ))] contains the same type of terms as 
$o (Q 2 ,Q 2 i), but each is a power of N higher. Thus [^(Q 2 , Q 2 ) — \n.(a s {Q 2 ) / a s {Q 2 ))} is subleading 
to <£,}((3 2 , Qj), and indeed &f(Q 2 ,Q 2 ) is a renormalization-scheme-dependent quantity, as it 
must be in order to absorb the change in &q(Q 2 ,Q 2 ) resulting from a change in the definition 
of a s (Q 2 ) under a change in renormalization scheme when working beyond leading order. Hence, 
we should factor exp[$^(Q 2 , Qj ) - \n(a s (Qj) / a s (Q 2 ))] out of the first term. Since [3>j~ (Q 2 , Qj) - 
ln(a s (Qj )/a s (Q 2 ))] is of the same form as the zeroth-order-in-7V part of &f(Q 2 ,Qj), then it 
should also be factored out of the leading-order expression for Fl(N,Q 2 ). Thus, we are left with 

'36 - 8N f ~ 



Fl(N,Q 2 ) =^Ml(F L ,o(N)+F L ,o(N,Qj) - 



+ 



2vr 

a s {Qj) 
2vr 



27 



F 2 , (iV))exp[$+(Q 2 ,Qf)] 



36 - 8N f 
27^ 



(4.54) 



F 2 , (N). 
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Using the one-loop running coupling constant, as is appropriate for a leading-order expression, the 
whole of (4.54) is not only manifestly factorization scheme independent, but also renormalization 
scheme independent, as one would hope. Also, Fl^(N,Q 2 ) is constructed precisely so as to make 
the expression unchanged, at this order, under a change in starting scale. Hence, within this 
expansion scheme (4.54) is genuinely the leading-order expression for Fl(N,Q 2 ). 

We now turn our attention to the more phenomenologically important case of the 
structure function F 2 (N,Q 2 ). For technical simplicity we begin with the expression for 
a 2 (Q 2 )(dF 2 (N,Q 2 )/da s (Q 2 )). Using the definitions introduced in our discussion for the longi- 
tudinal structure function we may write (4.44) as 

° U ° \F L AN)+ho(N,Ql)-( S -^^)F 2 , (N))>< 



2Nf \C 9 Ll (a s {Q 2 )/N) J ' uv ' "' uv V 27 

exp[<D+(Q 2 , Qj) + *+(Q?, Q 2 )] - a s {Q 2 )^^F 2fl {N) exp[<D r (Q 2 , Q])] 
+ higher order in a s and/or N. 

(4.55) 

The factorization scheme independence of this complete expression guarantees that the term 
(l} g (a s (Q 2 )/N) + Cf jl (« s (Q 2 )/iV)7° 9 (as(Q 2 )/iV))(Ci !li0 /Ci !l (a s (g 2 )/iV)) is a factorization- 
scheme-invariant quantity. Indeed, it was shown by Catani and Hautmann that ( ; yj g (a s (Q 2 )/N) + 
C^ l (a s (Q 2 )/N)^ g (a s (Q 2 )/N)) and C 9 L ^(a s (Q 2 )/N) could each always be expressed in terms of 
the product of a factorization-scheme- and renormalization-scheme-independent factor (which 
they calculated) and a scheme-dependent factor, where the scheme-dependent part is the same for 
both [18]. Using their results it is a trivial matter to find that 



( 7 } 9 MQ 2 )/A0 + C| jl (a s (g 2 )/iV) 7 (? (a s (Q 2 )/iV))). 



3 / C 9 



L,1,0 



2N f \Cl l {a s {Q 2 )/N), 
1/3 °° \ (456) 



2vrV2 

which is clearly both factorization scheme and renormalization scheme independent. Thus, our 
expression for the explicit part of (4.44) is entirely in terms of factorization-scheme-invariant, and 
hence physically meaningful quantities, i.e. 

~ a '^ 2 ^ ~^^^ =^ s {Q 2 )llL{^s{Q 2 )/N){F Lfi {N) + F Lfi {N,Q 2 I ) - (^^^jF 2j o(N))x 

exp[$+(Q 2 ,Q 2 ) + cI»+(Q 2 ,Q 2 )]- as (Q 2 )l^F 2! o(iV)exp[<I>r(Q 2 ,Q 2 )] 
+ higher order in a s and/or N. 

(4.57) 
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This expression is analogous to that for Fl(N,Q 2 ) in (4.53), except that it has still not been 
explicitly separated into inputs and evolution terms. In order to do this we must rewrite (4.57) as 

- a ^ Q2) d 2%V =^(QW2L(as(Qj)/N)(F L AN)+~F L AN,Qj) ~ (^^)*2,o(iV))x 

exp[$+(Q 2 , Qj) + *+(Ql Q 2 ) + $+(Q 2 , Q 2 )] 



4 2Nf -a s (Ql)F 2 , (N) exp [* r (Q 2 , Qj) - ln(||^|) 
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where &f(Q 2 , Qj ) is a series of the same form as <3?+((3 2 , Q 2 -), defined as 



+ higher order in a s and/or N, 

(4.58) 



^(Q 2 ,Q?) = - / i-nrlnKtfJitJV.a^? 2 ))^^? 2 ). (4.59) 

Ja s (Q 2 ) aa s {q ) 

It is now clear how we obtain the "leading part" of this expression, i.e. factor out the 
subleading parts of the evolution, [<frf(Q 2 , Qj) + 4>| (Q 2 Qj)] (we note that <&t(.Q 2 >Q 2 i) 1S a 
renormalization-scheme-independent contribution to this subleading evolution), and [&i(Q 2 , Qj) — 
ln(a s (Qj) /a s (Q 2 ))]. This leaves us with the leading-order expression 



dF 2 (N, Q 2 " 
da s (Q 2 , y 



(F L , (N) + F Lt0 (N, Qj) - ( 36 27 8iV/ ) F 2f0 (iV)) exp(d>+(Q 2 , Qj)) - a s{ Qj)i^LF 2fi {N). 
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(4.60) 

Again, this expression is insensitive to changes in starting scale, up to higher order, and using the 
one-loop coupling constant, is renormalization scheme independent as well as factorization scheme 
independent. 

Finally we consider the expression for F 2 (N,Q 2 ) itself. It is now a relatively simple matter 
to write this in terms of factorization-scheme-independent quantities and in terms of inputs and 
evolution terms. Using the definitions we have already made and defining &^(Q 2 , Qj) by 

^<n 2 n^ P (Q2f) d i ( a ^hM q 2 )/N) \ 

$i(Q ,Qi = -/ ~~j — r^T ln 5-7^7 — r^V\ — ) da s(q )■ 4.61) 

J a s (Q 2 ) da s (q 2 ) V T gg {N,a s (q 2 )) J 

we may factor terms into inputs and evolutions. Also, by construction it is guaranteed that the 
change of (4.43) at Qj under a change in starting scale will cancel the change in the evolution in 
the first term under a change in starting scale up to higher orders. It is also clear that (4.43) at 
Qj is both renormalization scheme and factorization scheme independent, as we require, as long 
as fi(N,Qj) is scheme independent. The requirement that if Qj = All the input reduces to the 
nonperturbative input F 2j0 (N) determines ff{N, Qj) uniquely. It must be the scheme-independent 
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quantity -a^illMQ]) /N) /f gg {a s {Q 2 ) /N)){F Lfi {N) - ( l^oW)- Hence, the in- 

put for F 2 {N,Q 2 ) is 



F 2 (iV,Q 2 ) = F 2 , (-/V) + Qs (Q 2 -) Zf^'^.^l F L , (N, Qj) + higher order in a s and/or TV, (4.62) 
and the expression for F 2 (N,Q 2 ) is 

F 2 (iV,Q 2 ) = a.(^) %^^g^ (^,o(jN0 +^,o(iV,Qf) - f ^^WiV^x 

eX P [c&+(Q 2 , Q 2 ) + <D+(Q 2 ) Q2 } + ^+ (Q 2 ; q 2)] + p 2 ^ N) exp[$r (Q2 > 

+ higher order in a s and/or N. 

Looking at (4.63) it is clear that there is only one formally LO (in our expansion scheme) 
input multiplied by a LO evolution, and that is F 2t $(N) multiplying unity. But this is obviously 
completely independent of a s , and is rather trivial. At next-to-leading order, or equivalently, at 
leading-a s -dependent order, we include the whole of (4.63) except that we factor exp[<3?^(g 2 , Q 2 ) + 
&t(Q 2 i Q 2 i)]i out of the first term. Hence, at leading-a s -dependent order we have 

F 2 ,o(N, Q 2 ) = F 2: o(N) exp[$- (Q 2 , Qj)] 

+ a ^^ ^|(t(S ) /% ) ( ^'° (iV) + ^' o(iV ' Qj) " F2AN)) exp[ ^ +(g2 ' 

n (0 2, ll L MQj)/N) ( p (m ( 36-8N f \ 

" as{Ql) i° g9 MQ 2 )/N) { Fl < o( N) -{-^r- ) F2AN)) - 

(4.64) 

In this expression there are clearly no terms which mix if we were to make a change in definition 
of the coupling a s — > a s + ea 2 , and hence we can consider it as a leading-order expression. If the 
one-loop coupling is used, it is both factorization scheme and renormalization scheme independent. 

We now have the full set of LO expressions in the leading-ln(l/x) expansion scheme. We 
could obtain the correct scheme-independent expressions for the structure functions at higher 
orders within this expansion scheme, but we choose to finish at leading order. The labour required 
to obtain higher-order expressions becomes progressively greater and we would obtain expressions 
requiring unknown anomalous dimensions and coefficient functions. Working to NLO we would 
need to calculate all the NLO evolutions and all the NLO inputs. This would require all the 7g a 's, 
7 2 a 's and C" 2 's. There is optimism that these NLO terms will soon be known [36], and once this 
is so the full NLO scheme-independent expressions should be calculated. 

We should make some comments about our LO scheme-independent expressions for Fl(N, Q 2 ), 
F 2 (N,Q 2 ) and (dF 2 (N, Q 2 )/dln(Q 2 )). First we note that all depend on factorization-scheme- 
independent combinations of the 7° a 's, 7} a 's and Cf^s (along with the input parton distributions, 
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where go(N,Q 2 ) is also factorization scheme dependent). There is, however, unlike in the loop 
expansion no terribly simple prescription for how one uses the anomalous dimensions and coefficient 
functions in order to arrive at the expressions. One must simply examine the form of the solutions 
and factor out subleading parts. 

By comparing (4.62) with (4.54) we now see that there is a very direct relationship between 
the inputs for our two structure functions at small x, i.e. we have a definite prediction, up to 
additive nonperturbative parts which are flat at small x, for one in terms of the other, as well 
as approximate predictions for the form of each. Also comparing with the form of (4.60), these 
two inputs for the structure functions are directly related to the slope of (dF 2 /dln Q 2 ) for small 
x at Qj. We do not yet know at what Qj it is most appropriate to choose the inputs, but it is a 
nontrivial requirement that the inputs for the three expressions are of the correct form and related 
in the above manner at any Q 2 . 

We should also make some mention of why factorization scheme dependence can be very large 
in this expansion scheme. In order to do this let us consider (4.42) as an example. A representative 
example of the way in which factorization-scheme-dependent calculations are done is to consider 
this expression evaluated with ^ g {a s {Q 2 )/N), ^f\{a s {Q 2 )/N) and C 9 L 1 (a s (Q 2 )/N) known in some 
particular factorization scheme, but ; y gg (a s (Q 2 )/N), which is unknown, either set equal to zero, or 
guessed by imposing some ansatz such as momentum conservation. As a first comment we consider 
the input. The input in terms of parton distributions is multiplied by C£ l {a s {Q 2 ) /N), a series 
which, in general, diverges at A = X(Q 2 ). Under a transformation of the type (4.29) this series will 
be multiplied by (U gg (a s (Q 2 )/N))~ 1 . A number of scheme transformations that are considered 
have the series U g (a s {Q 2 )/N) also becoming singular at N = X(Q 2 ) (e.g. [24][25]). This leads to 
powerlike growth of the form x _1_A< - < ^ / ^ as x — > 0, but the magnitude of the powerlike behaviour 
and the manner in which in which it is approached depends on the the precise behaviour of the 
coefficients in the series. Hence, changes in factorization scheme can lead to very marked differences 
in the form of the gluon at small x, or if the gluon is kept roughly constant (e.g. it is attempted 
to predict the form of the input structure function by assuming a form for the input gluon), to 
significant changes in the form of Fl(x,Q 2 ). Very similar considerations also hold for F2(x,Q 2 ) 
because the input depends strongly on Cf i (a s (Q 2 )/N), which transforms in the same way as 

c£,i,o(«.(Q?)/tf)- 

Examining (4.42) we can also see how the evolution may be strongly factorization scheme 
dependent. As we have already mentioned, the whole of &f(Q 2 ,Qj), including ~f gg (a s (Q 2 )/N), 
must be used in order to obtain a factorization-scheme-independent expression. The integrand in 
§i(Q 2 , Q 2 ) is a series which is a power of a s (Q 2 ) down on the series j gg (a s (Q 2 )/N). However, in 
many popular factorization schemes the coefficients in the incomplete, or incorrect series for this 
integrand are much larger than those in j gg (a s (Q 2 )/N) (helped by the fact that many of the early 
coefficients in ^f gg (a s (Q 2 ) /N) are zero), e.g. they commonly behave roughly like (12 In 2/7r) ra n -3 / 4 , 
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whereas the coefficients in ^® g {a s {Q 2 )/N) behave roughly like (12 In 2/ir) n n~ 3 / 2 . Hence, the in- 
correct $f(Q 2 ,Q 2 ) can have a dominant effect on the evolution. Under changes of factorization 
scheme the coefficients in the incorrect series can change by amounts similar to their own magni- 
tude. Therefore, the evolution of the structure functions in terms of a given input can appear to 
have a very strong factorization scheme dependence. Once again, this is true for the evolution of 
F 2 (x, Q 2 ) as well as for Fl(x, Q 2 ): the influence of the incorrect &f(Q 2 , Q 2 ) can be more important 
than that of ^(Q 2 , Qj) and "/ 2 l(N, a s (Q 2 )) combined, where for the latter the coefficients in the 
series are again relatively small. 

Alternative calculational procedures can lead to expressions which are not only similar to 
(4.42)-(4.44) with incorrect or missing 7* Ja s (Q 2 )/N), but which have additional factorization- 
scheme-dependent terms. These will be formally of higher order than the terms in (4.42)— (4.44), 
but again can have very large coefficients in the series expansions. This can lead to even more 
dramatic effects than those outlined above, e.g. [22]. 

Once jg g (N, a s (Q 2 )) is known in a given scheme and (4.42)-(4.44) can be calculated correctly 
there is no guarantee that the correct &f(Q 2 , Q 2 ) is not larger than ^(Q 2 , Q 2 ). If this is the case, 
^i(Q 2 ,Qj) will then have a large, but at least definite, effect. However, because it is a formally 
NLO correction to the structure functions, if ^(Q 2 , Q 2 ) is introduced then the full set of NLO 
expressions, both evolution and input factors, must be calculated at the same time: the correct 
calculational method respects renormalization scheme independence as well as factorization scheme 
dependence. This requires many more terms than just &f(Q 2 , Q 2 ). Hopefully, the complete NLO 
expression, as well as being factorization scheme independent, will also cause only fairly small 
changes to the LO expressions. 

Having obtained our full set of leading-order expressions, we can also examine the difference 
between F 2 (N,Q 2 ) and (dF 2 (N,Q 2 )/dlnQ 2 ) in this expansion scheme. There is a distinction 
between F 2 (N,Q 2 ) and (dF 2 /d InQ 2 ) even in the usual loop expansion. Differentiating a fixed- 
order expression for F 2 (N,Q 2 ) and using the /3-function evaluated to the appropriate order in a s 
results in the fixed-order expression for {dF 2 /d InQ 2 ) plus terms of higher order in a s (Q 2 ), which 
depend on the /3-function. The size of these extra terms is of the same order as the renormalization 
scheme uncertainty, and the distinction between the fixed-order expressions for F 2 (N,Q 2 ) and 
(dF 2 /d InQ 2 ) is of similar magnitude to the distinction between renormalization schemes. The 
distinction when using the small-x expansion appears more graphically, and is not only dependent 
on terms in the /3-function. If we differentiate (4.64) with respect to a s (Q 2 ) we obtain 



-«?(Q 2 ) 



dF 2j0 (N,Q 2 ) 
da s (Q 2 ) 



2 



=a s (Q 2 )j 1 2L (a s (Q 2 )/N)(^ 
(F L:0 (N)+F Li0 (N,Q 2 )- 




-n s (Q 2 )- 



4 2N 



F 2fi (N)exp[^(Q 2 ,Q 2 )}. 



9 67T&, 



'0 



(4.65) 
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This is clearly not exactly the same as (4.60), the difference being due to the additional 
terms {7 gg (a s (Q 2 )/N)/^ g (a s (Q 2 j )/N)) in the first term, and the factors (a s (Q 2 )/a s (Q 2 j)) and 
exp[$^(Q 2 , Q 2 )\ in the second term. All of these factors are unity at the boundary of the evo- 
lution, and the two expressions are therefore identical in this limit, i.e the inputs are the same. 
Therefore, it is the evolution terms which are different when comparing (4.60) and (4.65). Writing 



7 ° gg (a s (Q 2 )/N)\ r n iQj) d 



^ g (a s (Q 2 )/N)J ^ 



a s (Q ) = a s (Qj)exp 



(Q2 ) da s {q 2 )""'" 
a s (Qj 



HllMq 2 )/N))da s ( q 2 ) 



In 



(4.66) 



(4.67) 



MQ 2 ) 

we see that the terms present in (4.65) but absent in (4.60) are NLO evolution terms. Thus, as in the 
loop expansion, the difference between the fixed-order expression for (dF^/dhi Q 2 ) and the InQ 2 
derivative of the fixed-order F 2 consists of terms of higher order. However, in the leading-ln(l/x) 
expansion this difference exists between even LO expressions, and even at leading order we have 
to decide which of the two expressions to use (although there is very little difference between the 
choices in practice). It may be argued that in certain senses {dF 2 {N ,Q 2 ) / d\nQ 2 ) is more natural 
because it is a "real" perturbative quantity, beginning at first order in a s , as does Fl(N,Q 2 ). We 
will discuss and confirm this choice in §4.4. 

We should also make some comment about longitudinal momentum conservation. The first 
moment of the parton distributions is interpreted as the fraction of momentum carried by that 
type of parton, and it is usually required that f s (l,Q 2 ) + g(l,Q 2 ) = 1. For this to be true for 
all Q 2 then, of course, (d(f s (l,Q 2 ) + g(l, Q 2 ))/dhiQ 2 ) = 0, and from the renormalization group 
equations this is true if 

7// (l, a s (Q 2 )) + 7. 9 /(l, a s (Q 2 )) = 0, lfg (l, a s (Q 2 )) + lgg {l, a s (Q 2 )) = 0. (4.68) 

When expanding the anomalous dimensions order by order in a s (Q 2 ), it is easy to specify that (4.68) 
be true for the anomalous dimensions at each order, and to define a wide variety of factorization 
schemes which maintain this. This guarantees that the fraction of momentum carried by the nth~ 
order parton distributions is conserved at each n. Sometimes, all the momentum is designated to 
be carried by the zeroth order part of the solution, but this need not be the case. Most often the 
input is implicitly assumed to be the all-orders input. In this case it is true that it must carry all 
the momentum, but this method destroys the strict ordering of the solution. 

The situation is not as simple for the leading-ln(l/x) expansion. As can be seen from the 
matrix in (4.34), the leading-order 7 contains entirely positive entries for N = 1, and is clearly not 
consistent with momentum conservation: /q(Q 2 ) carries a constant amount of momentum while 
that carried by go(Q 2 ) is constantly increasing with Q 2 . In a general factorization scheme there is no 
reason that working to a finite higher order will restore the relationship (4.68). Two general methods 
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have been proposed to restore momentum conservation [21], [20], [37]. Both methods guarantee 
momentum conservation only if one truncates the series for the 7's at the 7 l5 s and solves the 
whole renormalization group equation using this truncated 7. Hence, both prescriptions destroy 
any sense of ordering the solution correctly. We advocate that the most sensible approach is to 
obtain a well-ordered solution for the structure functions, i.e. regard a correct treatment of the 
physical quantities as of paramount importance. Thus, we simply take the hint offered us by the 
zeroth-order anomalous dimension, and accept the fact that momentum is not conserved order by 
order in this method of expansion. How badly it appears to be violated, however, will depend very 
much on which factorization scheme we claim to work in, and we can choose the violation to be 
negligible. We will discuss this in more detail in §4.5. 

Finally, we can also discuss the relationship between our scheme-independent solutions and 
the ones which would be obtained using Catani's physical anomalous dimensions. One can solve for 
the LO structure functions using these effective anomalous dimensions in exactly the same way as 
we solved for the parton distributions in the previous subsection. In the same way that we have the 
relationships between the LO anomalous dimensions, (1.1) and (4.27), we have relationships between 
the effective anomalous dimensions as seen in (3.7). Using these relationships it is straightforward 
to obtain the analogous expressions to (4.36), (4.37) and (4.38): 

F 2 (N,Q*) =a s (Q>) ^^^M { F L , (N) + ~F L ,o(N, Qj) - (^^)^o(iV))x 



exp 



JaJQi) ««(r) a s {q 2 ) 



4 2JVj 



^ + higher order in a s and/or N, 

(4.69) 



F L (N,Q 2 ) ={h,o(N)+P L , {N,Q 2 i) ~ ( m J Nf Ko(iV))x 



exp 



27 

f^M^/AT) - ( ^1 )fi L (a s (g 2 )/iV) 



p^) f° L (a s (g 2 )/AQ + — ^ V 27 y— — daJq 2 ) 

L S (Q 2 ) a2 s(l 2 ) ' a s {q 2 ) 

4 2Nj 

+ ( 36 ~27 Nf ) F2 '° {N) {^W)) +highCT ° rder in ° s and/ ° r N ' 

(4.70) 
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and 

- al{Q2) d ^m =^(Q 2 n L MQ 2 )/N)(F L>0 (N)+~F L!0 (N,Qj) - (^^)F 2 , (N))x 



cxp 



/ TLL[a ° q 2 > /N) + da.tf) 

2Nf (a (0 2 )\ '^o 
- I a s (Q 2 )F 2 , (N)— f + higher order in a s and/or TV, 

(4.71) 



where F Lfi (N,Qj) is a function of rather than "fgg- 
We see that, once we make the identifications 



•99 

and 



7° (a s (Q 2 )/N) = r° LL (a s (Q 2 )/N), l\ L MQ 2 )/N) = T\ L {a s {Q 2 )/N), (4.72) 



= f i L MQ 2 )/A0 - (^|^)f k(« s (Q 2 )/iV), 



(4.73) 



(4.70) is identical to (4.53), (4.71) is identical to (4.57) and (4.69) is identical to (4.63). Of course, 
the identifications (4.72)-(4.73) are exactly what we obtain from the definitions of the physical 
anomalous dimensions in §3. 

Thus, we are able to reach these expressions for the structure functions somewhat more di- 
rectly by using the physical anomalous dimensions, and do not have to worry about problems with 
factorization scheme dependence (though we do have to calculate the physical anomalous dimen- 
sions in terms of known coefficient functions and anomalous dimensions of course). Once we have 
obtained these expressions using the physical anomalous dimensions we may then separate each 
of the terms into input parts and evolution parts (where more of this, but certainly not all, has 
already been done automatically when using the physical anomalous dimensions) and keep the 
most leading parts, obtaining once again the LO expressions (4.54), (4.60) and (4.64). So, using 
the physical anomalous dimensions leads us in a rather more direct manner to the correct leading- 
order expressions. If we were to work to higher orders, the amount of simplification obtained by 
using the physical anomalous dimensions rather than working in terms of parton densities would 
increase significantly. 

However, we stress that one will always automatically obtain factorization-scheme- 
independent answers by working to well-defined orders in physical quantities even when working in 
terms of partons. Also, we stress very strongly that even if one uses the physical anomalous dimen- 
sions, care is still needed to obtain expressions which are consistent with renormalization scheme 
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dependence, and there is no simple prescription for obtaining the correct results even in terms of the 
physical anomalous dimensions. If we were to solve the evolution equations for structure functions 
using Y° L L , r° 2 , Y\ L and T\ 2 we would obtain the explicit parts of (4.69)-(4.71) (with T 1 LL = 0) plus 
corrections of higher order. Terms must still be factored out of these expressions in order to obtain 
the true leading order structure functions. We will discuss the phenomenological consequences of 
this in §6. 

In this section, we have derived well-ordered, and consequently, factorization-scheme- 
independent expressions for structure functions in the leading-ln(l/x) expansion (which should 
be useful at small x) up to the order which is useful at present. This expansion does, however, 
sacrifice any attempt to describe the structure functions at large x. We would hope there is some 
expansion scheme which will be useful at all values of x. In the next section we will show that there 
is indeed an expansion scheme which satisfies this criterion, and argue that it is the only really 
correct expansion scheme. 

4.3. The Renormalization Scheme Consistent Expansion. 

In order to devise an expansion scheme which is useful at both large and small x we would a 
priori expect that we would need to use anomalous dimensions and coefficient functions at low orders 
in both a s and the leading-ln(l/x) expansion. There have already been various methods along these 
lines, however, they have all been scheme dependent. We demand complete consistency of our 
expressions for physical quantities with renormalization scheme invariance (which then guarantees 
factorization scheme independence). Consequently, our approach will be different from those used 
previously, and the results and conclusions will also be somewhat different. 

To begin, let us consider what we have meant by "consistency with renormalization scheme 
dependence" so far in this paper. In both the loop expansion and the leading-ln(l/x) expansion we 
demanded that once we had chosen a particular renormalization scheme and chosen to work to a 
particular order in this renormalization scheme then we would include all terms in our expressions 
for the structure function which were of greater magnitude than the uncertainty due to the freedom 
of choice of renormalization scheme (i.e. the uncertainty in the definition of the coupling constant), 
and no others. In each case the leading-order term consisted of the lowest-order inputs multiplying 
the lowest-order evolution terms. If working with the n-loop coupling constant the uncertainty in 
its definition is of order a™ +1 . Thus, the uncertainty of the input or evolution when working to 
nth-order is the change in the leading-order input or evolution if a s — > a s (l + ea"). Hence, the 
uncertainty in the whole structure function is of the order of the change of the leading-order part 
under such a change in the coupling. Therefore the n t h~order renormalization-scheme-independent 
expression includes all complete terms smaller than this change. 
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This definition does give us a well-defined way of building up an ordered solution to the struc- 
ture functions, but relies upon the definition of a given expansion scheme. It leaves an ambiguity 
about how we define the leading-order expressions and in how we define the order of terms com- 
pared to this leading-order term. Our two examples, the loop expansion, where the size of a term 
is determined simply by its order in a s , and the leading-ln(l/x) expansion, where ln(l/x) is put on 
an equal footing to a s , are just two examples of expansion schemes. Both have potential problems: 
in the former one does not worry about the large-ln(l/x) terms which can cause enhancement at 
small x of terms which are higher order in a s , and in the latter one does not worry about the fact 
that at large x, especially as Q 2 increases, it is the terms of lowest order in a s that dominant. 
Hence, both should have limited regions of validity. 

The shortcomings of these two expansion schemes come about because, even though any given 
order contains no terms which are inconsistent with working to the same given order in a particular 
renormalization scheme, in neither case does it include every one of the terms which are consistent 
to working to a given order in the renormalization scheme. In each expansion scheme some of the 
terms appearing at what we call higher orders are not really subleading in a s to any terms which 
have already appeared. Thus, despite the fact that for a given expansion method these terms are 
defined to be the same order as uncertainties due to the choice of renormalization scheme, they are 
not terms which can actually be generated by a change in renormalization scheme. 11 

In order to demonstrate this point more clearly we consider a simple toy model. Let us imagine 
some hypothetical physical quantity which can be expressed in the form 

oo oo oo oo / f f^ilW 3 

H(N,a s (Q 2 )) = a s (Q 2 ) £ a mn N n ee J>UQ 2 ) £ bJ ^jf 1 ) , (4.74) 

m=l n= — m i=0 j = l — i ^ ' 

where the expansion in powers of N about N = is convergent for all N. The first way of writing 
H(N,a s (Q 2 )) as a power series corresponds to the loop expansion, where we work order by order 
in m, out to m = k, and use the /c-loop coupling. The second corresponds to the leading-ln(l/x) 
expansion where we work order by order in i, out to i = I, and use the (I + l)-loop coupling. 
Let us, for a moment, consider the LO expression in the loop expansion, a s (Q 2 ) J2^L-i a inN n . 
The coupling is uncertain by 0(a 2 (Q 2 )) and hence the uncertainty of the leading-order expression 
(i.e. the change due to a change of the coupling) is ~ a 2 (Q 2 ) Yl^L-i bi n N n . There is no change 
with powers of ./V less than —1, and hence any such term is not really subleading. Similarly, the 
uncertainty of the leading-order expression in the leading ln(l/x) expansion contains no terms at 
first order in a s (or with positive powers of N), and such terms are not really subleading either. The 
full set of terms contained within the combination of both leading-order expressions is genuinely 
leading order, and is therefore renormalization scheme independent by definition. 

11 Similarly, they cannot be generated by a change in renormalization scale. 
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Perhaps the best way in which to write our expression for H(N, a s (Q 2 )) in order to appreciate 
these points is 

oo oo oo oo 

H(N, a s (Q 2 )) = Y, Nm E c ^a n s (Q 2 ) + £ N~ m £ c mn a n s (Q 2 ), (4.75) 

m= — 1 n=l m=2 n=m 

i.e. as an infinite number of power series in a s (Q 2 ), one for each power on N. Each of these 
series in a s (Q 2 ) is independent of the others, and the lowest order in a s (Q 2 ) of each is therefore 
renormalization scheme independent and part of the complete LO expression for H(N,a s (Q 2 )). 
The full LO expression for H(N,a s (Q 2 )) is therefore 

oo oo 

H (N,a s (Q 2 ))= Y, ^c^asiQ^ + YcmmN^aTiQ 2 ) 

(4.76) 

^ s (Q 2 )£a 0n iV« + £s{^Y. 

n=-l j=2 ^ ' 

Hence, the combined set of terms considered LO in the two expansion schemes comprise the full 
set of renormalization scheme invariant, and thus truly leading-order, terms. By considering 
H(N,a s (Q 2 )) written in the form (4.75), and considering how the coefficients in the expression 
must change in order to make the whole expression invariant under a redefinition of the coupling 
constant, a s (Q 2 ) — > a s (Q 2 ) + 0(a™(Q 2 )), we see that the n t h-order expression for H(N, a s (Q 2 )), 
which should be used with the n-loop coupling constant, consists of the sum of the first n terms 
in each of the power series in a s (Q 2 ). Thus, the full n t h _ order expression consists of the n t h^order 
expression in the loop expansion plus additional terms with inverse powers of N greater than n. 

Similar arguments have already been applied to the anomalous dimensions and coefficient 
functions, for example [18] and particularly [20], though in a somewhat weaker form. Here we 
take a strong, inflexible viewpoint and insist that the complete renormalization-scheme-consistent 
expressions (with no artificial suppression of leading— In (1/x) terms [20]), must be used. Further- 
more, and very importantly, the expressions used must be those for the physical structure functions, 
not for the factorization-scheme- and renormalization-scheme-dependent coefficient functions and 
anomalous dimensions. 

When considering the real structure functions the situation is technically a great deal more 
complicated than our toy model, but the principle is exactly the same. This can be seen by 
examining the the LO expressions for the structure functions in the two expansion schemes already 
considered. There is some overlap between the two, but each contains an infinite number of terms 
not present in the other. Considering the change of each leading-order expression under a change 
of coupling of 0{a 2 ) we see that, as with our toy model, the change in the LO structure functions 
in the loop expansion contains no terms in the LO expressions in the leading-ln(l/x) expansion, 
and vice versa i.e., none of the terms contained within each of the LO expressions are generated 
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by uncertainties at higher order in the other. Therefore, they should really all be regarded as 
genuinely LO, and be included in the full "leading order" expressions for the structure functions 
which use the one-loop coupling constant. Since these expressions contain all the parts of the one- 
loop expressions, and also contain leading-ln(l/x) terms as well, they should be able to describe 
the data over the full range of parameter space (except very low Q 2 , of course), as we would like 
from our correct LO expressions. However, there are two main complications when considering 
structure functions in comparison to our simple toy model. One is that the structure functions are 
combinations of evolution parts and input parts, rather than one simple power series in a s (Q 2 ). 
The other is that in general the physical anomalous dimensions, out of which the perturbative 
parts are constructed, are nonanalytic functions which cannot be expressed as power series about a 
particular value N for all N. They have singularities at N = 0, and also at negative powers of N 
(as well as possible a s (Q 2 )-dependent nonanalyticities due to resummation effects, e.g. the branch 
point in T° LL (N,a s (Q 2 )) at N = X(Q 2 )). We will deal with this second complication first. 

Let us consider the perturbative parts of the expressions for the structure functions. The 
singularities at negative integer values of ./V mean that we cannot write any physically meaningful 
quantity as just a power series about N = (or about (N + 1) for the nonsinglet case). Any such 
power series expansion will have a radius of convergence of unity, and a series expansion which 
applies over the whole range of N does not exist. In order to overcome this problem let us consider 
making the inverse transformation of some physically relevant perturbative function A(N, u s (Q 2 )) 
to x-space. The inverse of the Mellin transformation (2.1) is 

A(x,a s (Q 2 )) = — e W [CN]MN,a s (Q 2 ))dN. (4.77) 

where the line of integration is to the right of all nonanalyticities, and £ = ln(l/x). Since 
A(N, a s (Q 2 )) has, in general, singularities for all nonpositive integers, this whole integral may 
be evaluated by performing an infinite series of integrals, each with a contour centred on a given 
singularity, and not extending as far as unity from this singularity, i.e. not reaching any of the 
other singularities. Within each of these contours the function A(N, a s (Q 2 )) may be expanded as 
a power series about the singularity, i.e. we may write 

oo 

A(x,a s (Q 2 )) = — Y,x n exp[Z(N + n)]A n ((N + n),a s (Q 2 ))d(N + n), (4.78) 

n=0 J c ™ 

where A n ((N + n),a s (Q 2 )) denotes A(N,a s (Q 2 )) expanded as a power series about N = — n. The 
integrals will produce functions of £, which do not sum to integer powers of x, and each of the 
integrals in (4.78) will be independent and physically relevant in its own right. Hence, we must 
consider the complete moment-space expression as an infinite number of expressions of the form 
(4.75), each one having power series expansions in terms of (N + n), where n = — ► oo. The 
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expression for each n is then related to the part of the x-space expression behaving ~ x n . Of 
course, in practice, unless we want to examine the details of the structure function for x very close 
to 1, we can ignore all n greater than a finite, relatively small constant. 

Thus, when calculating the expressions for the perturbative part of the singlet structure func- 
tions, we are only concerned about LO terms beyond lowest order in a s for the specific case of 
n = 0. For n > we take the whole LO expression to be the one-loop expression. The terms 
we ignore by making this necessary decision are those which are LO in ln(l/x) at first order in x. 
Although these terms grow like a™ ln 2m_1 (l/x), there is no evidence that their coefficients are any 
larger than those for the zeroth-order-in-x logarithms. Since the resummed terms at zeroth order 
in x only begin to make an impact as x falls to ~ 0.1 (as we will see), and only become dominant 
for x much smaller than this, the effect of terms like ia™ \n 2m -\l/x) should be very small in 
comparison. Indeed, the effect of those terms of the form ia™ ln 2m_1 (l/x) which are actually 
known, i.e. m = 2, are indeed negligible. In a similar manner, we only consider the one-loop 
expressions for the nonsinglet structure functions in practice: the other LO parts of the expressions 
again lead to small-x enhancement of the form ia™ ln 2m_1 (l/x), which is very small compared 
to the singlet small-x enhancement, and there is only detailed data at very small x for the total 
structure function. 

We now return to our first problem, i.e. the fact that the structure functions are expressed in 
terms of both inputs and evolution parts. We discuss the case of the nonsinglet structure functions 
as an example first. 

4.4. The Renormalization Scheme Consistent (RSC) Nonsinglet Structure Functions. 

We consider a nonsinglet longitudinal structure function. For the nonsinglet structure func- 
tions the physical anomalous dimensions contain no singularities at N = 0, so the leading-ln(l/x) 
behaviour comes from singularities at N = — 1. Expanding about N = — 1, the full LO physical 
anomalous dimension can be written in the form, 



Tl s (N + l,a s m) =«,(Q 2 



M^ + i) 

m= — l 

m=l vv 1 / / m= i 



(4.79) 

a s {Q 2 , 



.(* + !) 

The first sum is just T^^N+l), the one-loop anomalous dimension expanded in powers of (N+l). 
The second sum contains the leading singularities in (N + 1) for all other orders in a s {Q 2 ). The 
final sum is included because, despite the obvious fact that it is a power of (N + 1) _1 down on 
the second sum, a series of this form cannot be created from the second sum by a change in the 
definition of the coupling of 0{a 2 s {Q 2 )). Therefore, the third sum is not subleading in a s (Q 2 ) to 
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the second sum, and must be renormalization scheme independent. Integrating (4.79) between Q 2 
and Q 2 , and including the overall power of a s (Q 2 ) for the longitudinal structure function we obtain 
the leading-order evolution 



[r (Q? YvV 1 h ± \ f ^LY\ da s (q 2 

exp il (Q2) \S\NTi bm + Cm ){WTW) ) 



(TV +1)2; Jb a 2 (q^ 



. (4.80) 



In order to construct the full LO input we must consider how the evolution term (4.80) changes 
under a change in starting scale, and therefore how the input must change in order to compensate 
for this. With a little work we can see that the full LORSC input is 

a s (Qj)F» s (N)(c»l l (N + 1) + C£f fJ (tf = -1) [\u{Q 2 / A N L s )a s {Q 2 )x 

(4.81) 

where 6o = a -i and Q) = ®- The ^ rs ^ term is just the lowest-order input in the loop expansion, 
while the second includes all the leading-ln(l/x) terms. 

Hence, the full LORSC expression for the nonsinglet longitudinal structure function, expanded 
about N = —1, is 

F? S rsc,o(N,Q 2 ) = « s (g?)F 2 7(iV)(c^(iV + 1) + C^N = -l)(ln(Ql/A» s )a s (Q 2 )x 

m=0 VV ' J / m=0 n=0 W 'J 



2 
7 

iV + 1 ^— ' " m V (iV + l) 2 J ' Z^V-™ 1 2 ^-n-m-n— ^ (jy + 1)2 



(iv+1) 2 ; ;& a 2 (g2)j- 

(4.82) 

In practice, we will only use the one-loop expression for the nonsinglet structure functions for all 
the singularities in the anomalous dimensions. This is because of the phenomenological reasons 
given at the end of the last subsection, and also because of lack of knowledge of the full physical 
anomalous dimensions. The terms ~ a™(N + l) 2 ™ -1 i n the anomalous dimension are all known 
[38], however, there is little knowledge yet of the terms of the sort ~ a™(N + l) 2r ™- 2 . These are 
an intrinsic part of the LORSC expression for the nonsinglet structure function, and should be 
calculated and included in order to give a true indication of the effect of leading-ln(l/x) terms. 
Hence, we believe that calculations of the nonsinglet [23] [39] (and polarized [40]) structure functions 
which claim to include leading-ln(l/x) corrections are presently incomplete. 

We now consider the nonsinglet structure function FJ? S (N,Q 2 ). This leads us back to our 
previous question of whether we should use the full RSC expression for F 2 (N,Q 2 ) or that for 
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(dF 2 (N, Q 2 )/d In Q 2 ). In order to illustrate the difference between the two we consider the simpler 

nonsinglet case, and also pretend for the moment that there is no small-x enhancement at higher 

/ 2 \r£ s i( Ar+i) 

orders in a s . Hence, the LO term in the evolution is just ( ^ \q' 2 ) ) ■ ^ ne hiput may be 

written as a power series in a s (Q 2 ), as we saw in §4.1, and is of the form 

F£o(N, Qj) = F 2 N S (N)[1 + a s (Qj ) ln(Q 2 /^ s )r^(iV + 1) + higher order in a s (Q 2 )]. (4.83) 

Therefore, as well as the lowest order a s -, and hence Q 2 -- dependent part of the input there is 
also the "sub-lowest-order", (^-independent part. These two terms are clearly of different order, 
but under a change in renormalization scheme, both remain unchanged and both should therefore 
appear in the LO definition of the structure function. This mixing of orders seems rather unsat- 
isfactory, and comes about because for F 2 NS (N,Q 2 ) the structure function still exists at zeroth 
order in a s and hence, it is not a perturbative quantity in quite the sense way as F^ S (N,Q 2 ) 
or (dF 2 NS (N,Q 2 )/dlnQ 2 ), both of which vanish with a s . So in our simplified model the LORSC 
expression for FJ? 3 (N,Q 2 ), obtained by combining the LO input and evolution, is 

FSiscoW' Q 2 ) = F 2 N o(m + HQj/A™)T^ l (N + 1)] (^gfj) ■ (4-84) 

This consists of two parts which are clearly of different magnitude, and it is clear that the same 
effect will be seen for the singlet structure function. 

If we instead consider (dF 2 NS (N,Q 2 )/d InQ 2 ) the full expression is 

dF2 Zl N Q2 Q2) =T 2 , N s(N + l,a s (Q 2 ))F 2 NS (N,Qj)x 

! Q (r 2 , NS (N +l,a s (q 2 ))- -A^\n(F 2iNS (N + l,a s (q 2 ))))d In q 2 . 

(4.85) 

Hence, the input may be written as 



exp 



^ (X ' Q1 ' ' F»s {N) [ as (Q 2 )r 2 '°' l (N + 1) + a 2 M)(HQVA?s )iT ^ {N + ^ 



d\nQ 2 ; q2 

+ r 2 '^(iV + l)) +0(a 3 s (Q 2 ))]. 

(4.86) 

The 0(a 2 (Q 2 )) piece is renormalization scheme dependent, so this time we have a LORSC input 
which is of a given order in a s (Q 2 ). The full LO expression for (dF 2 (N, Q 2 )/d\nQ 2 ) is then 

/dF, N n s (N 2 )\ ( n, (n2\\^2:li(N+i)-i 

( £ %^)r f - Wa - < « (yv + 1) (Sil) ' (487) 
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This is rather more satisfactory than the renormalization-scheme-independent expression for 
F 2 NS (N,Q 2 ) itself (4.84), and hence we choose (dF 2 (N, Q 2 )/d\nQ 2 ) to be the perturbative quantity 
we calculate, in both the nonsinglet and singlet case. 

If we wish to calculate the structure function F 2 (N, Q 2 ) itself to a given order we will do this 
by integrating the given order expression for (dF 2 (N, Q 2 )/d\nQ 2 ) between Q 2 and Q 2 , and adding 
it to F 2 (N,Qj) evaluated to the same order. For example, in our simplified nonsinglet model we 
would integrate (4.87) and add this to the explicitly written part of (4.83). This results in the 
effective LO expression 



+ a.(Q?) HQ 2 i/A^)T^(N + 1) 



(4.88) 



Of course, this whole discussion of F 2 NS (N,Q 2 ) has been rather simplified by the assumption 
that the higher-order-in-a s terms in the physical anomalous dimension do not contain higher 
singularities in (N + 1). Recognizing that they do, we obtain 



7V + 1 



m=l v v ' ' m=0 

exp 



i+e~ =1 m< +1 (q 2 )/(^+i) 



.i+e:=iLk +1 w?)/(,v + d 2 ' 

and 

PgiscfiW Q 2 i) = F 2 N o S W (l + ln(Q? /^ 5 )a s (Qf ) (V 2 %(7V + 1)+ 



(4.89) 



7V + 1 ^ "\(N + 1, , 

m=l ^ y ' ' m=U n=0 

(4.90) 

where for F 2 NS (N, Q 2 ) the coefficients in the series in (N + l) _1 are not necessarily the same as for 
F^ S (N,Q 2 ), hence the slightly different notation. However, a_i = a_i and b m = b m (probably 
[39]). The c m are not equal to the c m though, so there is no guarantee that F 2 N ^ SC (x, Q 2 ) and 
Fl,rsc,o( x i Q 2 ) wu ^ behave in the same way in the small-x limit. 

When all the b m and c m are known they can be used to present an argument for the form 
of the small-x behaviour of nonsinglet structure functions. Until this happens our discussion of 
the LORSC calculation of the nonsinglet structure functions is rather academic. However, it has 
enabled us to discuss many of the issues in a simpler framework than if we had gone directly to the 
singlet structure functions. We will discuss these singlet structure functions next. 
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4.5. The Renormalization Scheme Consistent Singlet Structure Functions. 

When calculating the singlet structure functions we cannot just construct the complete LO 
evolution and input and combine these to obtain the LO expression. Each of the component parts 
of the LO expressions for F2(N, Q 2 ) and Fl(N, Q 2 ) (we omit the superscript S in this section) must 
consist of LO input parts and evolution parts, but it is not obvious what these are. In order to find 
the full LORSC expressions for the singlet structure functions we will have to work in steps. We will 
consider only the full LO expression with the perturbative factors expanded about the particular 
value of N = 0, and the simplest way to proceed is to solve the evolution equations in terms of 
physical anomalous dimensions and structure functions. We have already proved that in the loop 
expansion the LO expressions only depend on the one-loop physical anomalous dimensions, and 
in the leading-ln(l/x) expansion the LO expressions depend only on (3.7). Hence, it is only the 
combination of these anomalous dimensions which is considered in our solution. 

We cannot simply write the physical anomalous dimension matrix 

a s (Q 2 )F° £ l L (N) + Tl L ( as (Q 2 )/N) a s {Q 2 )Tf 2 (N) + r° L2 (a s (Q 2 )/iV) , 



a s {Q 2 )Tf L {N)+a s {Q 2 )T\ L {a s (Q 2 )/N) a s (Q 2 )T^(N) + a s {Q 2 )T\ 2 (a s {Q 2 ) /N) . 

(where T° L L {a s {Q 2 )/N) = T° LL {a s {Q 2 ) /N) with the one-loop component subtracted out, etc.), 
and solve the renormalization group equations because the full solution contains terms which are 
not properly of leading order. We must choose some way of solving for the structure functions 
systematically which enables us to extract the true LO behaviour in as simple a manner as possible. 
In order to do this we take account of the fact that the one-loop solutions for Fl(N,Q ) and 
(dF 2 (N,Q 2 )/dlnQ 2 ) must be part of the complete LORSC solutions. Hence, we split our anomalous 
dimension matrix up into the form 

(Tf L (N)+ T%(N)\ ( Tl L (a s {Q 2 )/N) r° L2 (a s (Q 2 ) /N) \ 

«.W ) y T i{N) r 0,* (iy)y ! + { as{ Q2 )T l L{as{Q 2 )/N) as{ Q2 )T l 2{as{Q 2 )/N) ) > I J 

and solve by treating the second matrix as a perturbation to the first. Doing this we obtain the 
one-loop solutions as the lowest-order solutions and can systematically calculate corrections to 
this, extracting the parts of these "corrections" which are leading order. 

So, first let us consider the solution to the 0(a s ) renormalization group equation with boundary 
conditions F^ l {N,Q 2 j) = F L (N) and F 2 ' l (N,Q 2 T ) = F 2 (N). We may write the solution for the 
longitudinal structure function as 

/ rr>2\\ f o,i,+ (A0 / / n 2\\ f°- l -(N) 

F^N,Q 2 )=F^ + (N)[^M^ +^'"( iV )(^y) ■ ( 493 ) 

where r°' ; '+' - (N) are the two eigenvalues of the zeroth-order physical anomalous dimension matrix 
(which are the same as the eigenvalues of the zeroth-order parton anomalous dimension matrix), 
and in practice 

F°' l > + (N) = F L (N) - 36 ~ 2 * Nf F 2 (N) + 0{N), F° L ' l '-(N) = 36 ~* Nf F 2 (N) + O(N). (4.94) 
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Having chosen to write the lowest-order solution for the longitudinal structure function in this way 
we may then write the lowest-order solution for F 2 (N, Q 2 ) as 

F^\N,Q 2 ) = e^N)F^(N)^^j +e - {N) ^- {N) (^M"j , (4.95) 

where 

e+(N) =( r °'" +( ^ ) " jv r "' i(iV) ) = N/6 + 0(N 2 ), 

ni ( 496 ) 

e-(N) JI^M^m) = ( 27 ^ + (N) 
The first corrections to the one-loop solution may be obtained by solving the equations, 



d (H\N,Q 2 )\_ 2 ( F°£ l L (N) T%(N)\ (F?(N,Q 2 )\ 
dlnQ 2 \Ff(N,Q 2 )) ~ asW '{^(N) T° 22 l (N)J \Ff(N, Q 2 ) J 



+ 



Tl L (a s (Q 2 )/N) -^llTl L (a s (Q 2 )/N) \ (F° L ' l (N,Q 2 ) 

K a s (Q 2 )Tl L (a s (Q 2 )/N) -^EA as {Q 2 )T\ L (a s (Q 2 ) /N) ) \F°' 1 (N,Q 2 ) ) ' 

(4.97) 

where we have used the relationships between the physical anomalous dimensions in order to sim- 
plify the second matrix. We proceed as follows. First we define the vectors 

(4.98) 

and write 

F cl (N,Q 2 ) = e + (N)F cl ' + (N,Q 2 ) + er(N)F cl '-{N,Q 2 ). (4.99) 

We also define projection operators p + (N) and p~(N) (p + (N) ■ e + (N) = 1, p + (N) ■ e~{N) = 
etc.) which in practice are 

p+(N) = ( wj-m ) +Q(N), P~(N)=( 27 )+Q(N). (4.100) 

V 27 / V 8Af/-36 / 

Multiplying (4.97) by p + (N) or p~(N) leads to straightforward first-order differential equations 
for F cl '+(N,Q 2 ) or F cl ~{N,Q 2 ). Writing the zeroth-order solution F°' l (N,Q 2 ), in the form 

F°' l (N,Q 2 ) = e + {N)F°' l '+(N,Q 2 ) + e-{N)F°' l '-(N,Q 2 ), (4.101) 

these equation can now be solved by using the power series expansions of e +( ~~\N) and p +( -~\N) 
in terms of N. Only a small part of the overall solution contributes at LO. 
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Solving for the complete leading-order part of F + (N,Q 2 ) inductively in this manner, and 
making it formally insensitive to changes of Q 2 we get 



F L (N) ~ ( 36 27 8iV/ )F 2 (N))(eMHQ 2 i/A LL )Tl L (»s(Qj)/N)] - 1) 
exp 



a s (Q 2 ) 



\ If 

) exp U, 



'« S (Q 2 ) 



a 2 s {q 2 ) 



(4.102) 

We can solve for the corrections which are proportional to e~(7V) in exactly the same manner, 
using p~(N) instead of p + (N). Inserting the whole of F^ sc Q (N, Q 2 ) into the right-hand side 
the solution is relatively simple if we wish to keep only the LO parts. Adding to the solution at 
one-loop, the whole of the LO part of F~(N,Q 2 ) is 



? o,i,- 



(N) 



36 - 8N 



27 



27 



36 - 8N 



27 



a s (Q 2 )Tl L MQj)/N) - fTl L M Qj)/N)\\ ( a s {Q 2 
Tl L (a s (Q 2 )/N) 
i2\ /at\ _ N r b 



))\» S (Q 2 : 



"(AT) 



+ 



ti l mqiw) ) v l{n) - \-^-r {N) ) 

\ u sWJ/ Ua s (Q 2 ) oi 2 s {q 2 ) 

+ higher order, 



ex P [HQ 2 /A LL )Tl L MQj)/N)] (^0yj exp 



(4.103) 

where, in order to make the expression invariant under changes in starting scale, and also ensure 
that for Qj =A LL we have F~(N, Qj) = F®' l '~(N), we make the choice 



pcfull,- 



(N,Qj) 



36 - 8N f \ f a s (Q 2 )Tl L (as(Qj)/N) - f T LL MQ 2 )/N) 

rl L (a s (Q 2 )/N) 



27 



(f l (N) - ( 36 27 Nf ) F2{N) ) (eMHQ 2 i/A LL )Tl L (a s (Q 2 i)/N)} - 1). 

(4.104) 

We must now use the leading parts of F + (N,Q 2 ) and F~(N,Q 2 ) in order to obtain LORSC 
expressions for the structure functions. The way in which we have set up the calculation makes 
this very straightforward for the longitudinal structure function: we multiply F^ sc (N,Q 2 ) and 
Fft SC (N,Q 2 ) by (a s (Q 2 )/2ir), where the only part of F^ sc (N,Q 2 ) which contributes to the 
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LORSC expression for Fl(N,Q 2 ) is the one-loop part, 



"L,RSC,0 



(N,Q 2 



a s (Qj) 
2vr 



a s (Qj) 
« S (Q 2 ) 



(JV)-l 



exp 



fi L (a s (q 2 )/N) 



a s {q 2 



da s (q 2 ) 



F2' l ' + (N) + ( F L {N) 



36 - 8N f 



2 - jF 2 (N)J {eMHQ 2 i/A LL )Tl L (a s (Q])/N)] - 1) 



+ F°' ; ~{N) 



f°' l '-(N)-l 



(4.105) 

We also wish to find the LORSC expressions for (dF 2 (N, Q 2 )/d InQ 2 ) and for the input 
F 2 (N,Q 2 I ). We consider the former first. Using the form of e + (N) and e~(N) in (4.96) it is 
clear that, besides for the one-loop contributions, the LORSC expression will come from (27/(36 — 
87V/)) x (dF~(N,Q 2 )/d InQ 2 ) and from TV/6 x (dF+(N,Q 2 )/d InQ 2 ). Explicitly we obtain 



dF 2 (N,Q 2 

d InQ 2 J RSC0 



a s (Qj) 



e-{N)T > l '-(N)F£ l '-(N) 



0,1, 



a s (Qj 
a s (Q 2 ) 



2 f ' ! '-(Af)-l 



+ ( e +(A0r°><<+(A0Ff ' + (N) -r 2 ; L (iv)( f l (n) - ('^_|*> )f 2 (n) 



27 



+ r 1 2L (a s (Q 2 )/N)[F L (N) - ( 36 J Nf )F 2 (N) ) ex.p[ln(Q 2 j / A LL )T° LL (a s (Q 2 I ) /N)} 



27 



exp 



[/' 

J a 



T\ L {a s { q 2 )/N) 



(Q 2 ) 



a 2 (q 2 " 



da s {q 2 ) 



a s (Qj 
« S (Q 2 ) 



2 f ' ! ' + (Af)-l 



(4.106) 

where f (Tl L (a s (Q 2 )/N) - r| L (« s (Q 2 )/AT)) = a s (Q 2 ){T\ L (a s {Q 2 )/N) - v\ L (a s (Q 2 ) /N)) has 
been used and [N) ^F L (N) — ^ 36 ^ Nf ^ F 2 (N)^j is the part of the input which is common to 
both the one-loop and the leadmg-ln(l/x) input, and is subtracted from the former in order to 
avoid double counting. Similarly the input for F 2 (N, Q 2 ) is given by the relatively simple form 



F 2 , RS c,o(N,Qj) = F 2 (N) 

. rp 2 ^ L (a s (Q 2 )/iV) ( p 



(36 ~ 7 87V/) ^ 2 (iV)) (exp[ln(Q 2 M LL )r° LL (a s (Q 2 )/iV)] - 1) 
+ HQ 2 JA LL )a s {Q 2 )[e + {N)T^ + (N)Fl' 1 > + (N,Q 2 ) + e~ (N)T > 1 >- (N)F° L > 1 >- (N,Q 2 ) 

-r^ L (iV ) (F L (iV)-( 3 l^)F 2( iV))). 

(4.107) 

The third term is the renormalization-scheme-invariant order-a s (Q 2 ) input, which must compen- 
sate for changes in the one-loop evolution under changes in Q 2 . Again we explicitly extract a term 
oc ln(Q 2 /A)Y^° L (N) ^F L (N) - ( ^ 6 ~ 8 7 Nf ^j F 2 (N)^j in order to avoid double counting. The fact that 
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this term can be thought of as appearing from two different sources leads us to choose both our 
unknown scale constants in the input equal to the same value All- We note that our choice of 
inputs not only ensures Qf -invar iance up to higher orders, but are also such that our expressions for 
Fl,rsc,o(N,Q 2 j) and F 2i rsc,o{N,Q i ) reduce to the nonperturbative inputs if Q 2 = All- Having 
obtained our LORSC expressions for (dF 2 (N,Q 2 )/d InQ 2 ) and F 2 (N,Qj), then as already argued, 
in order to obtain our expression for F2(N,Q 2 ) we integrate (dF 2 (N,Q 2 )/d InQ )rsc,q from Q 2 
to Q 2 and add to the input F 2 ,Rsc,o(N,Qj). 

Thus, we have our complete leading-order renormalization-scheme-consistent expressions for 
the structure functions. These are significantly different from both the one-loop expressions and 
the leading-ln(l/x) expressions, although they do reduce to them in the appropriate limits. Indeed, 
once we include the 0(a s (Q 2 )) inputs for F 2 (N, Q 2 ) in the definition of the leading-order inputs in 
the loop-expansion (as we should), each of the terms in the inputs and evolution terms in (4.105), 
(4.106) and (4.107) contains a part which appears in both the LO expression in the loop expansion 
and the LO expression in the leading ln(l/x) expansion. Our full LORSC expressions are obtained 
rather more easily than above by simply letting each of the input and evolution terms become the 
combination of the terms in the two expansion schemes. In a sense this result is obvious, but it is 
necessary to verify this by deriving the expressions as above. 

Let us comment on the form of our final LORSC expressions. We note that as for the leading- 
ln(l/x) expansion, the LORSC expansion still leads to predictions for all the small-x inputs: 
predictions of each in terms of the nonperturbative inputs (which we imagine should be quite flat) 
and the nonperturbative scale All-, and also stronger predictions for the relationships between the 
inputs (although the scale Q 2 at which they should be imposed is not determined). 

Finally, we notice that each of the terms appearing in our expressions is manifestly renormaliza- 
tion scheme invariant, and it is clear that no terms are subleading in a s to any other terms, either in 
the input or in the evolution. If we had simply solved the renormalization group equations using the 
whole of the anomalous dimension matrix (4.91) then we would have obtained many terms which 
do not appear in our full leading-order expressions (4.105), (4.106) and (4.107). These would still 
be renormalization scheme independent, since our input anomalous dimensions are renormalization 
scheme independent. However, they would be of the same form as terms which are renormalization 
scheme dependent (e.g. we saw in 4.2 that the subleading-in-ln(l/x) evolution 3>+((3 2 , ) has 
a manifestly renormalization-scheme-independent part depending on T\ L {a s (Q 2 )/N) as well as a 
renormalization-scheme-dependent part depending on T LL (a s (Q 2 )/N)). These terms should be 
dropped, and (4.105), (4.106) and (4.107) are the correct expressions for the structure functions to 
be used with the one-loop coupling constant. 

46 



5. x Space Solutions. 

We shall now discuss how we use the expressions (4.105), (4.106) and (4.107) in order to obtain 
our expressions for the x-space structure functions and ultimately compare with data. The data on 
JF 2 (ar, Q 2 ) exist over a range of Q 2 from ~ 0.2GeV 2 to 5000GeV 2 , though in practice we will impose 
a lower cut on the data of Q 2 = 2GeV 2 , except for the HERA data where we choose Q 2 = 1.5GeV 2 
simply in order not to lose some of the very low x data. The threshold of heavy quark production is 
W 2 = Am 2 H , where W is the invariant mass of the hadronic system created from the struck proton 
(or neutron) and in the limit of zero proton (and/or neutron) mass is given by W 2 = Q 2 {x~ 1 — 1). 
tuh is the mass of the heavy quark. Thus, we clearly work in the limit where the up, down and 
strange quarks are effectively massless. However, we cross the 6-quark threshold of W 2 ps 20GeV 2 , 
and at the lower end of our range are in the region of the c-quark threshold of W 2 « 2GeV 2 . 
Hence, we will need expressions for the structure functions which cross quark thresholds. 

The correct treatment of these heavy quark thresholds is less well established than the treat- 
ment of effectively massless quarks, and is certainly more complicated. Here we use the prescription 
for treating heavy quarks outlined in [41]. This involves treating all the quarks as massless, but 
only allowing the heavy quarks to become active above the simple threshold Q 2 = m 2 H . Hence, the 
value of Nf appearing explicitly in any expressions changes discontinuously at this threshold. The 
running coupling constant, is defined to be continuous at the thresholds. It is determined by the 
relationship 



This complete prescription for treatment of the heavy quarks is consistent with the decoupling 
theorem, as it is guaranteed to provide the correct expressions far above or below any threshold: 
the increase in the coupling below a threshold compensating for the absence of virtual heavy 
quarks in calculations below this threshold. It is clearly rather unsatisfactory in the region of the 
threshold, with heavy quark structure functions having an abrupt threshold in Q 2 rather than a 
smooth threshold in W 2 . Work to rectify this is in progress, and will certainly involve the use of 
the heavy quark coefficient functions at leading order in ln(l/x) already calculated [13] [31]. 

We can now discuss the form of the x-space solutions for the structure functions. It is the 
perturbative part of the moment-space structure functions for which we can produce well-ordered, 
RSC expressions and the nonperturbative parts of the expressions for these structure functions, 




(5.1) 



where the central a s (Q 2 ) with Nf = 4 is defined by 



12vr 



(5.2) 



(33-2.(N f = 4))HQ 2 /A 2 QCDA Y 
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F 2 (N) and Fl(N), will be nonanalytic, complicated functions of N. Our complete moment-space 
expression for a general structure function (or derivative of a structure function) will be 



F i {N,Q 2 ) = P it2 



(N,a s )F 2 (N)+P itL 



(N,a 8 )F L (N), 



(5.3) 



where the P^ ( N, a s ) 's are the calculable perturbative components of the complete expressions 
which can be expanded as power series in both a s and N, i.e. examples of the physically relevant 
perturbative functions we discussed in §4.3. We obtain the x-space structure functions as follows. 
We factor out the perturbative part of each of these expressions which is proportional to either 
Fl(N) or F 2 (N), and take the inverse Mellin transformation of this perturbative function by inte- 
grating around a contour encircling N = 0, obtaining the leading-in-ln(l/x), at lowest-order-in-x 
parts of the perturbative P(x, a a )'s as well as all the lowest-order-in-a s parts at lowest-order- 
in-x. In order to obtain the leading-order-in-ln(l/x)-and-a s , at lowest-order-in-x, part of the 
perturbatively calculated structure functions we then convolute with the whole of the nonpertur- 
bative Fl(x) and F 2 (x) which are obtained by a complete inverse Mellin transformation of Fl(N) 
and F 2 (N). In principle, the full LORSC x-space structure functions are calculated by repeating 
this procedure for the LORSC moment-space structure functions where the perturbative parts are 
expanded about N equal to all negative integers, but in practice one only really needs include the 
one-loop, or leading-in-a s part of the full LORSC perturbative parts for all negative integers. 

In practice the calculation of the structure functions is performed in a rather different manner. 
The structure functions are calculated using a modification of the computer program that is used 
by Martin, Roberts and Stirling in their global fits to structure function data. This works in 
terms of parton densities directly in x-space. Input parton densities are specified at some scale 
Qo and the Q 2 evolution is calculated on a grid in x and Q 2 . This evolution is obtained by 
integrating up the renormalization group equations involving the complete parton distributions and 
full specified splitting functions. The structure functions are calculated by numerically performing 
the convolutions of the resulting Q 2 -dependent parton distributions with coefficient functions. 
Thus, in fact, we obtain the scheme-independent LORSC structure functions by working in terms 
of parton densities and choosing coefficient functions and splitting functions (i.e. an effective 
factorization scheme) which will reproduce the correct expressions as closely as possible. 

We will briefly describe our choice of effective splitting functions as follows. At first order in 
a s we choose the normal parton splitting functions. We then add corrections to these in order to 
reproduce our desired results. Denoting these corrections by A'y a i,(N, a s (Q 2 )), for a given parton 
to parton splitting function, and expressing these in moment space and in terms of previously 
discussed quantities for simplicity we get for jF L (x,(5 2 ) 



A lfg (N,a s (Q 2 )) = A lff (N,a s (Q 2 )) = 
A lgg (N,a s (Q 2 )) = \A lgf {N,a s {Q 2 )) = Tl L (a s (Q 2 )/N). 



(5.4) 
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This form is not too difficult to understand by looking at e.g. (4.105): the only leading-ln(l/a;) 
enhancement of the one-loop expression comes from the evolution in the longitudinal sector which 
is directly related to the enhanced gluon evolution and hence corrections to the gluon anomalous 
dimensions. For Tiix^Q 2 ) the corrections are a little more involved: 



A lfg (N,a s (Q 2 )) = a s {Q 2 )C 9 2Xl {N)v\ L (a s {Q 2 )/N) 
A Xff (N,a s (Q 2 )) = a s (Q 



V 27 3/ N f 2 > lA ' 



r\ L {a s {Q 2 )/N) / s 

(5.5) 



A lgg {N,a s {Q 2 )) = Tl L {a s {Q 2 )/N) - |A 7/s (7V, a s (Q 2 )) 
A lgf (N,a s (Q 2 )) = lTl L (a s (Q 2 )/N) - |A 7// (iV, a s (Q 2 )). 

The corrections to the quark anomalous dimensions are functions of Q 2 rather than Q 2 , ex- 
cept for the single power of a s (Q 2 ), because T\ L (a s /N) appears only in the input terms in 
(4.106). In the small-x, or small-iV limit we have the simplifications that A^f g (N, a s {Q 2 )) — > 
a s (Q 2 )(2N f /3)rl L (a s (Q 2 )/N) and A 7// (7V, a s {Q 2 )) § A 7fg (N, a s (Q 2 )). Thus, in this limit 
these corrections take on the standard form of a s (Q 2 )^j g (a s (Q 2 )/N) and a s (Q 2 )^jj:(a s (Q 2 )/N) 
minus their one-loop components. The prefactors multiplying these terms, which depend on the 
one-loop longitudinal coefficient functions, ensure that the leading-ln(l/x) enhancement of the 
rate of growth of J~2(x, Q 2 ) is directly coupled to the longitudinal structure, not to the gluon. This 
delays the enhancement to slightly smaller values of x, and reduces it a little. The contributions to 
the gluon anomalous dimensions beyond LO in ln(l/x) are present to counter the exponentiation 
of the corrections to the quark anomalous dimensions. 

We perform checks to see if the particular choice above does indeed lead to an accurate repre- 
sentation of the correct expressions. In order to do this we use the types of techniques outlined in 
[22] to transform the exact expressions (4.105), (4.106) and (4.107) to x-space obtaining analytic 
solutions for the structure functions. These expressions are very complicated, but are almost exact 
for very low values of x, and are compared, with the calculations performed using the computer 
program in x-space. For x ^ 0.1 our choice of the anomalous dimensions leads to agreement with 
the analytic expressions to an accuracy much better than the errors on the data in any appropri- 
ate range of parameter space. At higher x setting to zero all terms in the MRS program other 
than those coming from one loop, the expressions for the structure functions for x > 0.1, are very 
close to the one-loop expressions alone, which agrees with calculations taking the inverse Mellin 
transformation of the LORSC moment-space expressions. Therefore, the correct choice of effective 
splitting functions leads to the MRS program producing a very accurate approximation to our 
correct LORSC expressions for the structure functions over the complete range of parameter space. 

Now that it has been determined that our choice of splitting functions is correct, one can vary 
the input parameters in the standard way to obtain the best fit. In practice the inputs are of the 
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standard MRS form (4.24), in terms of partons, but these are constrained by demanding that the 
resulting structure functions at Q 2 are compatible with forms in (4.105), (4.106) and (4.107) (for 
some All) where the nonperturbative inputs are flat at small x, i.e. they must be described well 
by a function of the form 

Ti{x) = F i (l-x)^(l + e i x a5 + 7l x). (5.6) 

We note that with our choice of splitting functions and definitions of parton densities then 
momentum is not conserved by the evolution: in the best fit to ^{x^Q 2 ) discussed in the next 
section the total momentum carried by the partons at Q 2 = 2GeV 2 is 87% and at Q 2 = 5000GeV 2 
it is about 94%. Hence, the amount of momentum violation is at the level of a few percent. We 
have already defended this violation of momentum conservation in subsection 4.3. We now also 
point out that starting with our definition of the partons we could define new parton densities 
and splitting functions by defining non-zero C£ { ?\a s /N) and C f L { f{a s /N) beyond one-loop, and 
use the transformation rules in §2 to keep the structure functions unchanged. If these coefficient 
functions had negative coefficients then, compared to our prescription above, the low Q 2 parton 
distributions would need to be larger, and would hence carry more momentum. As Q 2 increased 
the effect of these new coefficient functions would decrease, and the extra amount of momentum 
carried by the new parton distributions compared to the original ones would decrease. Thus, the 
effect of such a redefinition of parton distributions would be to increase the amount of momentum 
carried by the partons at low Q 2 and also to slow the growth of momentum with Q 2 (or even to 
turn it into a fall). With a judicious redefinition of coefficient functions of this sort (with perhaps 
some dependence on a s (Q 2 ) as well as a s (Q 2 )) it should clearly be possible to find an effective 
factorization scheme where the momentum violation will be extremely small, and one may think of 
this as a "physical scheme" . We have not seriously investigated this redefinition of parton densities 
in any quantitative manner since it will not affect any physical quantities. 

Having specified precisely how we will perform our calculations of structure functions we can 
make some comment on the general form the structure functions have to take. The expressions 
for Tl(x,Q 2 ) and ^{x^Q 2 ) depend on T° LL (a s /N) and T\ L (a s /N). Both of these series have 
coefficients which behave like n _3 / 2 (121n2)/7r for very large n. This leads to a cut in the N- 
plane at N = a s 41n2 in both cases, and to the structure functions having asymptotic behaviour 
T{x,Q 2 ) ~ (lnx) _3 / 2 x _41n2as as x — > 0, where a s = a s (Q 2 ) and Q 2 > Q 2 . However, it was 
convincingly demonstrated in [20] (and discussed from a different point of view in [22]) that one 
only need keep a finite number of terms in the leading-ln(l/x) series if working at finite ln(l/x). 
In practice, if one works down to x ~ 10~ 5 and the nonperturbative inputs behave roughly like 
(1 — x) 5 , then keeping terms up to 10th order in the series for T LL (a s /N) and T\ L {a s /N) is more 
than sufficient. The series up to this order have the explicit form 

T° L (x) = x + 2.40x 4 + 2.07x 6 + 17.3x 7 + 2.01x 8 + 39.8x 9 + 168.5x 10 + • • • (5.7) 
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and 



2tiT£ l (x) = 1 + 2.5x + x 2 + x 3 + 7.01x 4 + 5.81x 5 + 13.4x 6 + 58.1x 7 + 64.7x 8 + 196.8x 9 + 650x 10 + • • • . 

(5.8) 

The low-order terms in both these series have coefficients which generally grow far less quickly 
than the asymptotic relationship a n+ i = 41n2a n ; some fall, and in the case of T° LL (a s /N) are 
even zero. Thus, in the range of parameter space we are considering we will have rather less steep 
behaviour than the asymptotic limit of JF(x,Q 2 ) ~ (lnx) _3 / 2 x~ 41n2Qs . Using the techniques in 
[22] one can derive analytic expressions for the structure functions in the small-x limit. Choosing 
Qj = 25GeV 2 and making a guess that A LL = 0.8GeV 2 we obtain ln(Q|/^4 LL ) m 3.4. We also 
choose nonperturbative inputs which behave like Fi(l — x) 5 for i = L, 2. As argued in [22], for 
values of x & 0.01 it is a good approximation to replace the (1 — x) 5 behaviour with the behaviour 
0(0.1 — x). 12 Doing this we can derive simple expressions for the form of the inputs for x 5s 0.01 
which will be accurate up to a few percent error. Using the series (5.7) and (5.8) and the expressions 
(4.105) and (4.107) we can take the inverse Mellin transformations, obtaining 



Jl(x,Q 2 ) «a s (Q 2 )/2vr 



F L + —F 2 )[1 + 3.4a s (Q 2 )£ + 5.8 v + 6.6 



and 



27 J V ' sv ^ /s 2! 3! 

+ 13.7< g y + 31.5^y + 56^fil! + 136^y) - ±F* 

(5.9) 

Q?) « 3Aa s (Qj)/2^F L + ±F 2 } (l + 4.2a s (Qj)C + 7.1 ( " s( g )g)2 + g.l^M^ 

+ 19 . 7 (M|M)! + 441 + M (*y + 206 (^M)O!) + F2 , 

(5.10) 

where £ = ln(0.1/x), and a s (Q 2 ) ss 0.2 if Kqcd,a rs lOOMeV. Putting in values of F L = 2.5 and 
F 2 = 1, choices which, as with the (1 — x) 5 behaviour, are roughly compatible with the high x 
data, we have a rough estimate of the form of the input structure functions at Qj = 25GeV 2 . 
The coefficients in the series in a s (Qj)^ multiplying (Fl + 4/27F 2 ) grow a little more quickly 
for ^(x^Qf) than for J^(x,Q 2 ). However, in the former case this contribution which rises with 
falling x is accompanied by the flat F 2 , whereas in the latter it is accompanied only by the nearly 
insignificant — 4a s (Q 2 I )/(54Tr)F 2 . Thus, for x « 10~ 4 Tl(x, 25) behaves approximately like O.lx -0 - 3 
but !F 2 {x, 25) is slightly flatter, behaving approximately like 0.6x~ 28 . These powers of x are clearly 
somewhat less steep than the asymptotic x -0 ' 5 . Comparison of the estimate of T 2 {x, 25) with the 



12 The argument of the 9-function depends on the power of (1 — x). Higher powers would require the 
step to occur at lower x and vice versa. 
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data in [1] and [2] shows a very reasonable qualitative agreement. Of course, there is no data for 
25) for values of x anything like this low. Being rather more general, we find that for any Q 2 
between lOGeV 2 and lOOGeV 2 , and with any sensible choice of A LL (e.g 0.2GeV 2 > A LL > 2GeV 2 ), 
then both F 2 (x,Q 2 j) and F L (x, Q 2 ) behave roughly like x~ - 3 for 0.01 > x > 0.00001. Hence, this 
type of behaviour can be taken as a prediction of the theory. 

One can also use the same techniques to make an estimate of the structure functions at values 
of Q 2 away from Q 2 . With the particular inputs above they do lead to good qualitative agreement 
with the data below x = 0.01, and we have every reason to feel encouraged by our results. However, 
the real test of our approach will be a complete global fit to the available data for Ti{x, Q 2 ) using 
the rather more accurate calculations, especially at large x, of the MRS program. We will therefore 
discuss these detailed fits next. 

6. Fits to the Data and Predictions. 

Once we have the general LORSC expression then by combining the singlet and nonsinglet 
components and varying all the free parameters (All, the soft inputs for (x, Q 2 ) and T% (x, Q 2 ) 
and the soft nonsinglet inputs), we obtain the best fit for the available JF 2 structure function data 
using a particular starting scale Qj. Choosing the renormalization scale to be Q 2 , the one-loop 
value for A^r /= 4 is fixed at lOOMeV, thus giving q s (M|) = 0.115. 13 This precise value is not 
determined by a best fit, but a value near to this is certainly favoured. 

There are some further details we should mention. Firstly, when obtaining a fit using the 
above approach, the values of Ti[x,Q 2 ') used are not precisely those published in [1] and [2]. This 
is because it is not Ti{x, Q 2 ) that is measured directly at HERA, but the differential cross-section. 
This is related to the structure functions as follows 



where y = Q 2 /xs « Q 2 /90000x at HERA and R = F L {x, Q 2 )/(J r 2 (x, Q 2 ) - F L {x, Q 2 )). Over most 
of the range of parameter space y is very small, Tl{x, Q 2 ) is likely to be small, and the measurement 
is essentially directly that of JF 2 (x, Q 2 ). However, at high y the value offoix, Q 2 ) must be extracted 
using some prescription for Tl{x,Q 2 ). Both HI and ZEUS, roughly speaking, obtain their values 
of Tl{x,Q 2 ) from predictions coming from NLO calculations of JF 2 (x,Q 2 ). Since the approach in 
this paper leads in practice to a somewhat lower prediction of Fl{x,Q 2 ), the values of JF 2 (x, Q 2 ) 
used in the fit must be altered to take account of this. Thus, the ^(xjQ 2 ) values are a little (at 
most about 6%) lower for the largest values of y than in those in [1] and [2]. 

13 Of course, since we are using a genuinely leading-order expression, any change in renormalization 
scale is exactly countered by a change in Ajv /= 4. However, it is encouraging that making the simple choice 
A*r = Q 2 leads to a value of a s (M%) which is nicely compatible with the usual value. 




(6.1) 
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The value of m 2 is chosen to be equal to 4GeV 2 in order to obtain a reasonable description of 
the data on the charm structure function coming from EMC [42] and from measurements at HERA 
[43]. The quality of the fit is shown in fig. 1. It is of a fair quality, with the predicted JT|(x,Q 2 ) 
perhaps being a little large in general at large x, and a little small at small x, a result which is 
qualitatively consistent with the fact that we have used a threshold at Q 2 = m 2 , rather than at 
W 2 = Q 2 (x~ 1 — 1) = 4m 2 . Of course, m 2 = 4GeV 2 is a little high compared with the values 
obtained from reliable determinations. We will comment on this value later. The strange quark is 
treated as being massless in our calculations, which is presumably a good approximation for the 
values of Q 2 considered. Also, we insist that the strange contribution to the structure function is 0.2 
of the singlet structure function minus the valence contributions (i.e. the "sea structure function" ) 
at Q 2 = m 2 . This ensures compatibility with the data on neutrino-induced deep inelastic di-muon 
production obtained by the CCFR collaboration [44]. 

Finally we consider the form of the gluon at large x. Within our effective factorization scheme 
we would expect the gluon distribution to be quite similar to that in the MS scheme at NLO-in-a<. 
for very large x. Thus, we demand that our gluon is qualitatively similar to that obtained from the 
WA70 prompt photon data [45] at x > 0.3, i.e. roughly of the form 2.5(1 -x) 6 at Q 2 = 20GeV 2 for 
values of x this large. Encouragingly, this is the type of large x behaviour that the best fit chooses 
for the gluon, and no strong constraint is needed. 

The fit is performed for a wide variety of data: the HI [1] and Zeus [2] data on ^^{x^Q 2 ) 
with 0.000032 > x > 0.32 and 1.5GeV 2 > Q 2 > 5000GeV 2 ; The BCDMS data [46] on ^ p (x,Q 2 ) 
with 0.07 > x > 0.75 and 7.5GeV 2 > Q 2 > 230GeV 2 ; the new NMC data [47] on ^ p (x,Q 2 ) 
and f% d (x,Q 2 ) with 0.008 > x > 0.5 and 2.5GeV 2 > Q 2 > 65GeV 2 ; NMC data on the ratio of 
F£ n (x,Q 2 ) to (x,Q 2 ) [48] with 0.015 > x > 0.7 and 5.5GeV 2 > Q 2 > 160GeV 2 , CCFR data 
[49] on F% N (x,Q 2 ) and ^ N (x,Q 2 ) with 0.125 > x > 0.65 and 5GeV 2 > Q 2 > 501.2GeV 2 ; and 
the E665 [50] data on f£ p (x,Q 2 ) with 0.0037 > x > 0.387 and 2.05GeV 2 > 64.3GeV 2 . One can 
see that the full range of data used in the fit covers an extremely wide range of both x and Q 2 and 
thus provides a very stringent test of any approach used to describe it. 14 

The result of the best fit to this data using the LORSC expressions with Q 2 chosen to be equal 
to 40GeV 2 , m 2 = 4GeV 2 and ml = 20GeV 2 is shown in table 1. Also the result of the fit to the 
small— a; data is shown in fig. 2. As one can see there is a very good quality fit to the whole selection 
of data, and thus over the whole range of x and Q 2 . Overall, the fit gives a x 2 of 1105 for 1099 data 
points. The fit shown is for the particular starting scale Q 2 = 40GeV 2 , but the quality of the fit is 
extremely insensitive to changes in this scale, as we expect from the method of construction of the 
solutions. The fit is essentially unchanged over the range 20 — 80GeV 2 , and we choose 40GeV 2 as 

14 We have not included ZEUS data for Q 2 > 2000GeV 2 . 
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the geometric mean. When Q 2 drops below 20GeV 2 the fit immediately gets markedly worse, due 
to the bottom quark threshold. 

There are 18 free parameters used in the fit: &-qcd,a (which we choose to describe as a 
parameter since, although it is fixed at lOOMeV, and small variations would no doubt improve 
the fit slightly, it can certainly not vary too much); four parameters for each of the valence quark 
contributions, which are described by functions of the form (4.24) where the normalization is set by 
the number of valence quarks; four parameters for the nonperturbative inputs for the two singlet 
structure functions, which are of the form (5.6); and the unknown scale All, where we allow All 
to be a free parameter for each Qj . We do not consider Q 2 as a free parameter, since it can take a 
very wide range of values. The parameter All, which we have argued should be a scale typical of 
soft physics, turns out to be 0.55GeV 2 for the fit starting at Qj = 40GeV 2 . This decreases a little 
as Qj increases and vice versa. For Qj = 40GeV 2 the soft inputs for the fit are roughly 

F[ (x) « 3(1 - x) 5 (l + O.lx - 5 - 0.2x), F£(x) « (1 - x) 3A (l - 0.65a; - 5 + 4.5a;). (6.2) 

These nonperturbative inputs lead to complete inputs of F[(x,Qj) « 3x -0 ' 33 and F 2 5 (x,Q 2 ) « 
0.65x-°' 3 for 0.01 > x > 0.0001, with the effective A increasing in both cases for even smaller x. 
Instead of forcing the nonperturbative inputs to be flat as x — > we could allow an asymptotic 
behaviour x~ x , where A Ss 0.08. This leads to an equally good fit. 

Thus, the LORSC fit seems to be a success. However, in order to gauge its true quality it is 
helpful to have some points of comparison, and we will discuss some alternative fits. The first we 
mention is that obtained by simply solving the evolution equations for structure functions using 
the full leading order physical anomalous dimensions (4.91), rather than using the more careful 
procedure described above. In this case if we drop the restriction on the relationship between the 
inputs and let them each have arbitrary small x form (the arguments relating the inputs do not 
really follow if we simply solve using the full physical anomalous dimensions) then the fit is nearly 
as good as the LORSC fit for Q 2 > 4GeV 2 . However, it gets markedly worse for Q 2 ~ 3GeV 2 and 
is extremely poor below this. If we do try to impose the relationship between the inputs (i.e. retain 
any explanatory power for the form of these inputs) , then the fit at small x is much worse than the 
LORSC fit, the value of (d^ix, Q 2 )/d \nQ 2 ) being rather too steep. Thus, phenomenologically the 
naive use of the LO physical anomalous dimensions is clearly inferior to the full LORSC calculation. 

We also compare to more conventional approaches. As in [30], the most recent MRS fits Ri 
and R 2 are shown. These are obtained using the standard two-loop method, where Ri allows 
A^£T 4 to be free (giving A^£r 4 = 241MeV) and R 2 fixes A^T 4 = 344MeV to force a better fit 

MS v MS ' MS 

to the HERA data. The MRS fits are useful for purposes of comparison for a number of reasons: 
the treatment of the errors in this paper is identical to that in the MRS fits, so consistency in 
this respect is guaranteed, any systematic differences due to differences in computer programs is 
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guaranteed to be absent, and the cuts in Q 2 for each data set are chosen to be the same in this 
paper as for the MRS fits. The number of free parameters in the NLO-in-a s fit is the same as 
in the LORSC fit. There is one more parameter in the NLO-in-a s fit due to the powerlike forms 
of the input gluon and singlet quark at small x being independent 15 whereas the small-x shape 
of the inputs for both structure functions is completely determined by the one parameter All va 
the LORSC fit. However, there is one less parameter in the NLO fit due to the normalization of 
the gluon being determined by momentum conservation. Essentially, the LORSC fit is a little less 
constrained at large x than the NLO fit, but rather more predetermined at small x. 

Comparing the LORSC fit and the MRS fits, it is clear that the LORSC scheme-independent 
fit is much better for the HERA data (even when compared to R 2 ), much better for the BCDMS 
data (even when compared to Ri) and similar in standard for the rest of the data. The overall fit 
is ~ 200 better for the whole data set; clearly a lot better. However, this direct comparison with 
the MRS fits is rather unfair: the MRS fits use the SLAC data [51] on T^{x^Q 2 ) which is not 
included in the LORSC fit due to much greater sensitivity of this data to potential higher twist 
effects than any of the other data sets, and also the fact that the quark and gluon distributions 
are input at lGeV 2 in the MRS fits is not helpful to their quality. If the best fit is obtained for 
the HERA data using the NLO calculation with massless quarks with inputs at Q 2 « 4GeV 2 , as 
in [27], and evolution performed downwards in Q 2 then the gluon becomes negative at very small 
x before Q 2 = lGeV 2 . Hence, starting at larger Q 2 where the gluon distribution is happy to be 
positive everywhere, produces a better fit. Thus, in order to make a more meaningful comparison 
of the LORSC fit with a NLO-in-a s fit we have ourselves performed a NLO-in-a s fit, called NLOi, 
allowing the normalizations to vary in the same way as in the LORSC fit, with exactly the same 
treatment of quark thresholds as the LORSC fit, with the input parton distributions chosen at 
Q 2 = m 2 and with the fit to exactly the same data. 16 As in the MRS fits the values of ^(x, Q 2 ) 
at small x are those quoted in [2] and [1] . The value of m 2 is chosen to provide a good description 
of the data on the charm structure function. The value needed is m 2 = 2.75GeV 2 , and the fit to 
the charm structure function data is shown in fig. 3. 17 The value of A— T 4 determined by the fit 

& MS 

is 299MeV. The quality of the fit is shown in table 1. 

15 We claim that the allowed powerlike form of the inputs at small x is against the spirit of a well-ordered 
perturbative expansion, as already discussed in §4.1. However, enforcing this rule would mean that the 
quality of the fit to the small-x data was very poor. Thus, we allow the more usual, unjustified choice 
for the small-x form of the inputs. In fact, at Q\ = lGeV 2 the gluon is strongly valence-like, while the 
sea-quark distribution ~ x~ ~ ' at small x. 

16 Once again we allow the small— x form of the inputs to be unjustified powerlike behaviours. At 
Q] — 2.75GeV 2 the gluon is quite flat while the sea-quark distribution ~ a; -1-0 ' 22 . 

17 As in the LORSC fit the choice of m 2 which gives a good description of the charm data is also the 
choice which gives the best global fit. 
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The higher starting scale for evolution has produced a better fit to the HERA data than the 
MRS fits. However, the best fit comes from allowing the normalization of the HI data to be at the 
lower limit allowed by the error on the normalization 18 , and is still not as good as the LORSC fit 
to the HERA data. Not including the SLAC data leads to a much better fit to the BCDMS data 
than the MRS fits, but again this is clearly not as good as the LORSC fit. The fit to the NMC data 
is much the same for the LORSC fit as for the NLOi fit, and the NLOi fit is a little better for the 
CCFR data. The overall quality of the NLOi fit is 1169 for the 1099 points, and thus is 64 worse 
than the LORSC fit. Therefore, the LORSC fit is still clearly better than the NLOi fit, but not 
nearly as convincingly as appeared to be the case when compared to the MRS fits. Nevertheless, 
it is encouraging that, while the overall fits to the relatively high x data, i.e. the BCDMS, NMC, 
E665 and CCFR data are similar in the NLOi fit and the LORSC fit, it is the fit to the small-x 
HERA data that is definitely better for the LORSC fit, as we would expect. This can be seen even 
more clearly if we examine the quality of the fit by separating the data into two sets: one where 
x < 0.1 and one where x > 0.1. This is shown in table 2, and demonstrates that the LORSC fit is 
superior at small x, while not quite as good as the NLOi fit at large x. 

This qualitative result is exactly what we would expect. The importance of the leading-ln(l/x) 
terms in the LORSC calculation can be quantitatively judged by how they affect the fit. If, after 
obtaining the best LORSC fit, all terms other than those in the one-loop expressions are set to 
zero, the quality of the fit is unchanged above x = 0.3, begins to alter slightly below this, and is 
clearly worse by the time we reach x = 0.1. Thus, the leading~ln(l/x) terms are important by 
this value of x. However, most of this effect is due to the terms at 0{a 2 s ), so the NLO expression 
at fixed order in a s should be insensitive to higher-order-in-a s leading-ln(l/x) terms down to x 
somewhat lower than 0.1. Therefore, above x = 0.1 the NLO fit should in principle be better than 
the LORSC fit since it contains terms at NLO in a s which are important at large x. However, 
the NLO fit should be considerably worse at small x since it does not contain many important 
leading-ln(l/x) terms. This is qualitatively in agreement with the comparison of the NLOi fit and 
the LORSC fit. 

However, the above comparison is somewhat incorrect because in the process of obtaining 
the best fit for all the data the NLOi fit may choose some parameters, particularly A^£r 4 , such 
that they mimic the effects of the leading-ln(l/x) terms, and a decent fit for the small x data is 
obtained to the detriment of the fit to the large x data. In order to check this hypothesis we have 
also performed a NLO-in-a s fit with A^L fixed at 250MeV, which we will denote by NL02- The 
results of this fit in terms of the different data sets is shown in table 1. The fit clearly improves 
compared to the NLOi fit for the BCDMS and CCFR data, and gives the best overall fit for the 
high x data sets. It worsens somewhat for the NMC data, and also for the HERA data, and overall 

18 The fit to the HI data continues to improve slightly for a normalization going down to 0.96. 
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is slightly worse than the NLOi fit, having a \ 2 of 1184 for the 1099 data points. Nevertheless, it is 
perhaps a truer representation of a real NLO-in-a s fit than NLOi since it gives a better fit to high 
x data, but not such a good fit to the small-x data which presumably require the leading-ln(l/x) 
terms. Whether one believes this argument or not, it is certainly clear that the LORSC fit does 
provide a better description of the data than any standard NLO-in-a s fit. 

This leads us to the question of the determination of a s (M§) using global fits to structure 
function data. The complete RSC expression for structure functions only exists at leading order, 
and this leaves the renormalization scale undetermined. Hence, a determination of a s (M§) does not 
really take place. It is not yet possible to extend the RSC calculation beyond the leading order due 
to lack of knowledge of NLO-in-ln(l/x) terms, but hopefully these will shortly become available 
[36], and when they do the NLO versions of (4.105)-(4.107) can be derived and put to use, and 
the NLO coupling constant determined. Until these full renormalization-scheme-consistent NLO 
expressions become available, we believe that it is incorrect to use the present NLO-in-a s fits 
to small-x structure function data in order to determine the NLO coupling constant. However, 
as already argued, the fixed-order-in-a*, expressions should be accurate for CCFR, BCDMS and 
NMC data (except perhaps at the lowest x values), which after all are still much more precise than 
HERA data, and fits to these data alone will provide the best determination of the NLO a s (Q 2 ). 

We have given what we hope are convincing arguments for our advocated approach for cal- 
culating structure functions. Not only do we claim theoretical correctness, and limited predictive 
power, but we also have good quality, very comprehensive fits to data on ^{x^Q 2 ). However, we 
are well aware that only further experimental tests can prove us right (or wrong). Hence, we now 
discuss our predictions for Tl{x, Q 2 )- 

So far we have only probed J-~l( x , Q 2 ) indirectly, i.e. it is simply related to the InQ 2 derivative 
of T2(x,Q 2 ) (as well as to the input ^(xjQ 2 )). Having tied down the nonperturbative inputs 
and All and Q 2 from our fit to J-2(x,Q 2 ), we have a prediction for J-l(x, Q 2 )- The result of this 
prediction for the fit with Q 2 = 40GeV 2 is shown in fig. 4, where it is compared to the prediction 
using the NLO-in-a s approach, in particular the NLOi fit- As one can see, it is smaller than the 
NLOi Fl{x,Q 2 ), but becomes steeper at very small x. The NLO2 fit gives a very similar form 
of J~l( x ,Q 2 ) to the NLOi fit- The LORSC prediction for !Fl{x^Q 2 ) is weakly dependent on the 
value of Q 2 chosen: the value at Q 2 = 5GeV 2 and x = 10 -4 varies by ±10% within our range of 
Q 2 (increasing with Q 2 ), and by less than this for higher x and Q 2 . 

Hence, measurements of Tl{x,Q 2 ) at x < 10~ 2 would be a good discriminant between fixed- 
order-in-a s calculations and those involving leading-ln(l/x) terms. However, the "determination" 
of Tl{x, Q 2 ) already performed by HI [52] is really only a consistency check for a particular NLO- 
in-as fit, and is by no means a true measurement of Tl{x, Q 2 ). In essence, all it proves is that the 
measurements of the cross-section are consistent with a particular NLO-in-a s fit to Fi(x, Q 2 ) when 
the relationship between the cross-section and the value of ^(xjQ 2 ) is determined assuming the 
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correctness of the NLO-in-a s expressions for both structure functions. This is a perfectly correct 
procedure, and should be adopted for any fit to JF 2 (x, Q 2 ) data, as it has been for the LORSC fit, 
but in itself says nothing about the validity of a different approach, or about the actual value of 
Fl(x,Q 2 )- Hence, real measurements of Tl{x,Q 2 ) are needed in order to find the real values of 
Tl{x,Q 2 ). From fig. 4 it is clear that such measurements at HERA would be an important (and 
probably essential) way of determining the validity of the approach in this paper, and the genuine 
importance of leading-ln(l/x) terms in structure functions. 

Another potentially important discriminant between different methods of calculating structure 
functions is the measurement of the charm structure function. Both the LORSC fit and the NLOi 
fit can provide good fits to the currently available charm data, as already seen, and the value 
of m 2 c in the NLOi fit is rather more satisfactory than that in the LORSC fit. However, these 
calculations involve incorrect treatments of the charm quark threshold. A more correct treatment 
of this threshold at NLO in a s [53] shows that the value of m 2 must be considerably reduced in 
order to produce the same sort of value of the charm structure function as the approach used in 
this paper. Indeed, as seen in [54], NLO-in-a s calculations with m 2 = 2.25GeV 2 undershoot the 
small-x data. A correct treatment of the charm quark threshold within the framework advocated 
in this paper has not been completely worked out (although, as already mentioned, work is in 
progress), and has certainly not been tested. However, a comparison of the theoretical calculations 
with the ever-improving data on the charm structure function seems potentially to be a very useful 
way of discriminating between different methods of calculation. 

In principle there are many other quantities which could be calculated within the LORSC 
framework and compared to those calculated using the NLO-in-a<. approach (or any other method) 
and to experimental data. Particularly obvious examples are the distribution of the transverse 
energy flow in the final state in lepton-hadron scattering and the the cross-section for forward jet 
production, for both of which there exists some experimental data which does not seem to be terribly 
well described by the order-by-order-in-a s approach, e.g. [55] and [56]. We have not performed 
a LORSC calculation of any such quantities, and it is very difficult to estimate the results, other 
than guess that they will probably lie somewhere between the fixed-order-in-a s predictions and 
those obtained using BFKL physics naively. Such calculations are obviously a priority, and work 
will begin soon. Only by comparing our theoretical predictions with a wide variety of experimental 
data can we determine unambiguously which theoretical approach is correct. 

7. Conclusion and Summary. 

In this paper we have derived expressions for the structure functions ^(x, Q 2 ) and Tl{x, Q 2 ) 
in a theoretically correct, well-ordered manner, within the framework of the renormalization group 
and collinear factorization. We have first done this for the particular expansions schemes which 
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order the expressions strictly in orders of a s , i.e. the standard loop expansion, or in terms of the 
leading powers of ln(l/x) for given power of a s , i.e. the leading-ln(l/x) expansion. In both cases 
we have demonstrated that a correct calculation in terms of structure function inputs and evolution 
automatically leads to factorization-scheme-invariant results which may be expressed in terms of 
physical anomalous dimensions. Thus, these physical anomalous dimensions are fundamental pieces 
in the correct expressions for the structure functions. However, we also demonstrate that in the 
case of the leading-ln(l/x) expansion the correct expressions are more difficult to obtain than in 
the loop expansion, and in both cases a correct calculation requires rather more care than just the 
use of these physical anomalous dimensions. 

We have then argued that both the above expansion schemes are restrictive, and lead to 
only part of the correct solution at any given order. We have shown that the only calculational 
method which is truly consistent with working to a given order in a s within a given renormalization 
scheme is the renormalization-scheme-consistent expansion, in which one works to a given order 
in both a s , and in ln(l/x) for given power of a s , for physical quantities, i.e. the inputs and 
evolution of the structure functions. Doing this one obtains the conventional, order-by-order- 
in-a s renormalization group results in the high x limit, and smoothly incorporates the leading 
ln(l/x) corrections. We have then derived the leading-order, renormalization-scheme-consistent 
(and hence renormalization scheme invariant) expressions for the structure functions Ti{x, Q 2 ) and 
^l(x,Q 2 )- In this derivation we have made use of the physical anomalous dimensions (which is 
not necessary, but is extremely convenient), but needed to do more work than simply include them 
in their leading form and solve the evolution equations. Indeed, this naive method of inclusion 
leads to a comparison with data which is far worse than that obtained from the correct LORSC 
calculation. 

As part of our overall approach we have also taken an unusual view of the starting scale Q 2 . 
Rather than trying to guess some form for the input at some particular Q 2 , or simply allowing 
the inputs at some arbitrary Q 2 to take any form they wish within a given parameterization, we 
have demanded firstly that our expressions should be formally insensitive to the choice of the input 
scale, and secondly that any deviation from the flat Regge-type behaviour of the structure functions 
must come from perturbative effects. Thus, our inputs take the form of nonperturbative functions, 
which are flat at small x, convoluted with functions of the physical anomalous dimensions which 
are evaluated at Q 2 , and determined by the requirement of insensitivity to the value of Q 2 . This 
leads to our inputs being determined entirely in terms of the flat nonperturbative inputs and one 
arbitrary scale All which roughly indicates the scale where perturbative physics should break 
down, i.e. A LL ^ lGeV 2 . This gives us a great deal more predictive power than more usual 
approaches. We have some idea of the form of individual structure functions at small x: for a range 
of sensible choices of Qj (lOGeV 2 - lOOGeV 2 ) and A LL we obtain f 2 (x, Q] ) roughly oc £-(0.25-0.33) 
for 0.01 > x > 0.00001, which is clearly in good qualitative agreement with the data. However, we 
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also have a much stronger prediction for the relationship between the small-x inputs. So as well as 
some predictive (or at the very least, explanatory) power for individual structure functions, once 
we choose the input for one, in practice T2(x,Q 2 ), we have determined, up to a small amount of 
freedom, the small-x inputs for d ^(x, Q 2 ) / dhiQ 2 and Tl{x,Q 2 ). 

Not only are the features of the LORSC calculation compelling, but they also work rather well 
in practice (for a range of Qj from 20GeV 2 - 80GeV 2 ). The LORSC expressions, including this 
constraint on the small-x inputs, lead to very good fits to the data, the \ 2 f° r the LORSC fit to 
1099 data on F 2 (x, Q 2 ) ranging from from 0.75 > x > 0.000032 and 1.5GeV 2 < Q 2 < 5000GeV 2 is 
better by at least 60 than any NLO-in-a s fit, a very significant improvement. This is even though 
the NLO-in-a s fit is allowed arbitrary, unjustified powerlike behaviour at small x, and the small-x 
inputs for Fi(x, Q 2 ) and dT 2 (x, Q 2 )/dlnQ 2 are largely independent. In fact, all of this superiority 
comes from the fit to the data with x < 0.1. Thus, we find the naively expected, but much refuted, 
result that the leading ln(l/x) terms are both important and helpful in fits to structure functions 
when included properly. As a caveat, it is certainly true that the calculations in this paper must 
be improved to take account of massive quark thresholds in a better manner (as is also the case 
with NLO-in-a s fits), and work towards this end is in progress. Nevertheless, with the present 
treatment we feel that the quality of the fit and the degree of explanatory (if not predictive) power, 
not to mention the theoretical correctness, give strong justification for the LORSC expressions. 

We do, however, recognize that the quality of the fit alone does not necessarily convince one 
that this approach has to be correct. In order to obtain verification we must exploit the factorization 
theory fully and compare with more and different experimental data. Hence, we have presented a 
LORSC prediction for Tl(x,Q 2 ), comparing it to that obtained using the NLO-in-a s approach. 
Hopefully there will be true measurements of J~l(x,Q 2 ) at HERA some time in the future with 
which we can compare these predictions. We stress the importance of such measurements to the 
understanding of the physics which really underlies hadron interactions. In the near future we will 
also have predictions of the charm contribution to the structure function within the framework 
of a correct treatment of the massive quark, and comparison with the ever-improving data on 
the charm structure function should also be a good discriminant between different theoretical 
approaches. Calculations and measurements of other, less inclusive quantities, such as forward jets, 
are another clear goal. 

Finally, our calculation is at present only at leading order due to the lack of knowledge of 
next-to-leading~order~in-ln(l/x) coefficient functions and anomalous dimensions, or equivalently 
of physical anomalous dimensions. This means that the NLO-in-a s approach is still in principle 
superior to our approach at high x, where the leading-ln(l/;c) terms at third order in a s and beyond 
are not important but the 0(a 2 ) terms are. Indeed, in practice the NLO-in-a s fits, are slightly 
better than the LORSC fit for data at x > 0.1, and we might expect any predictions coming from 
the NLO-in-a s approach to be more accurate than those coming from the LORSC approach down 
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to values of x somewhere in the region of 0.05. The lack of the NLORSC expressions also means 
that a true determination of a a (M\) from a global fit to structure function data is not yet possible, 
but that the best determination from fits to structure function data should at present come from 
using the NLO-in-a s approach, using only large x data. For a really fair comparison between the 
renormalization-scheme-consistent method and the conventional order-by-order-in-a s approach 
we really need the full NLORSC calculation. Hopefully the required NLO in ln(l/x) quantities will 
soon become available, and with some work we will be able to make such a comparison. We would 
expect that when using the NLORSC expressions the fit to the large x data would become at least 
as good as for the NLO-in-a s approach, and that the fit to the small-x data would be of at least 
the same quality as for the LORSC fit. We wait expectantly to discover if this is indeed the case. 
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Table 1 

Comparison of quality of fits using the full leading-order (including leading-ln(l/x) terms) 
renormalization-scheme-consistent expression, LORSC, and the two-loop fits MRSRi, MRSR2, 
NLOi and NL0 2 . For the LORSC fit the HI data chooses a normalization of 1.00, the ZEUS data 
of 1.015, and the BCDMS data of 0.975. The CCFR data is fixed at a normalization of 0.95, and 
the rest is fixed at 1.00. Similarly, for the NLOi fit the HI data is fixed at a normalization of 0.985, 
the ZEUS chooses a normalization of 0.99, and the BCDMS data of 0.975. Again the CCFR data 
is fixed at a normalization of 0.95, and the rest fixed at 1.00. Also, for the NLO2 fit the HI data 
is fixed at a normalization of 0.985, the ZEUS chooses a normalization of 0.985, and the BCDMS 
data of 0.97. Again the CCFR data is fixed at a normalization of 0.95, and the rest fixed at 1.00. 
In the Ri and R2 fits the BCDMS data has a fixed normalization of 0.98, the CCFR data of 0.95 
and the rest of 1.00. 
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Table 2 

Comparison of quality of fits using the full leading-order (including leading— In (1/a;) terms) 
renormalization-scheme-consistent expression, LORSC, and the two-loop fits NLOi an d NL02- 
The fits are identical to above, but the data are presented in terms of whether x is less than 0.1 or 
not. 
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Figure Captions 

Fig. 1. The description of the EMC and preliminary HI data for T% (x, Q 2 ) using the LORSC fit. 

Fig. 2. The curves correspond to the value of the proton structure function ^(xjQ 2 ) obtained 
from the leading-order, renormalization-scheme-consistent (LORSC) calculation at 12 
values of x appropriate for the most recent HERA data. For clarity of display we add 
0.5(12— i) to the value of Ti(x, Q 2 ) each time the value of x is decreased, where i = 1 — > 12. 
The data are assigned to the x value which is closest to the experimental x bin (for more 
details see the similar figure displaying the two-loop fits in [26] ) . E665 data are also shown 
on the curves with the five largest x values. The HI and ZEUS data are normalized by 
1.00 and 1.015 respectively in order to produce the best fit. 

Fig. 3. The description of the EMC and preliminary HI data for T^ix, Q 2 ) using the NLOi fit. 

Fig. 4. Comparison of predictions for JF L (x, Q 2 ) using the full leading-order, renormalization- 
scheme-consistent (LORSC) fit and the NLOi fit- For both sets of curves Tl{x,Q 2 ) 
increases with increasing Q 2 at the lowest x values. 
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